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ABSTRACT 


The first part of this thesis is principally a review of a 
number of geometrical formulations of the classical theory of rigid 
bodies. Chapter II describes a Lie group model of rotational motion 
originally due to Arnold [1], and Chapter III incorporates this model 
into the inhomogeneous symplectic mechanics corresponding to a general 
Lagrangian and state space. In Chapter IV these ideas are generalized 
to a homogeneous space-time formulation which includes the relativistic 
rigid test body of Kunzle [29] as well as a new description of a non- 


relativistic rigid body based on Galilean relativity. 


The second part of the thesis begins with a general discussion 
on constrained motion and then presents a geometrical theory of 
constraints. This is applied to a new formulation of constrained 
motion of a rigid body in curved space-time and several examples are 


presented to illustrate the methods. 
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CHAPTER I 


Introduction 


The description of rigid bodies and spinning particles has 
enjoyed a long and rich history ever since its inception with Euler in 
the earlier part of the 18th century. Euler's formulation was origin- 
ally given in the context of Newtonian mechanics; since then, particu- 
larly from the end of the 19th century to the present, these systems 
have been cast into the framework of the various new theories (and 
formalisms thereof) as they arose. To properly give a historical sur- 
vey would require the space of a book (indeed, see Corben [8] for a 
review of spinning particles and gyroscopes in special relativity and 


quantum mechanics) and is beyond the scope of the present work. 


We confine our historical remarks to noting: the emergence 
of a consistent set of equations of motion for a spinning test particle 
in a gravitational field with Mathisson [62] in 1937 and later 
Papapetrou [38] in 1950, the development of a multipole formalism for 
extended test bodies by Dixon [12], and the existence and uniqueness 
theorems for a centre of mass of an extended body (Beiglbock [5], 
Dixon [13], Madore [33]) leading to the equations of motion. All of 
these and similar derivations proceeded by integrating the conserva- 
tion law for the stress-energy tensor of an extended body. More 
recent approaches, such as Souriau [49], Bailey and Israel [3] 
(charged spinning particles in electromagnetic and gravitational 


fields), and Hanson and Regge [17] (special relativistic spherical 
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rigid body), are based on a Lagrangian and the associated variational 
principle. On a more geometrical level are the canonical symplectic 
formalisms: the canonical Lagrangian approaches of Kitinzle ({27], [29]) 
based on a presymplectic potential 1-form and the canonical Hamiltonian 


approach of Sniatycki and Tulczyjew [47] based on a symplectic 2-form. 


The purpose of this thesis is twofold: to survey the clas- 
sical formulation of the motion of a rigid body from a geometrical point 
of view and to develop a geometrical formulation of constrained rigid 


body motion in both Newtonian mechanics and general relativity. 


The contempory interest in rigid test bodies in general 
relativity seems to be due to, among other things, the gyroscope-satel- 
lite experiment that Schiff [43] proposed in 1960 to verify one feature 
of general relativity that hitherto had been unmeasurable (for example, 
see [6] and [14] for a discussion of the actual experiment, and [4], 
[35] and [36] for the theoretical background). This experiment, being 
carried out by C.W.F. Everitt and W.M. Fairbank and currently near 
completion, attempts to very precisely measure the precession Giatie 
spin axis of a rapidly spinning rigid test body. The various analyses 
on which the experiment is based are carried out, at least implicitly, 
in the context of the PPN ("parameterized post-Newtonian") formalism 
(see [35] for definitions and discussions of this formalism). In this 
framework, the effect of the earth's rotation upon the gravitational 
field shows up as off-diagonal terms in the PPN expansion of the metric 
tensor. It is only in the low-temperature, (almost) force-free envir- 
onment of the satellite experiment that one hopes to detect this extre- 


mely small effect. 
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Chapter II has the dual function of establishing the notation 
that will be used throughout as well as presenting a model of rotational 
motion that will be a key element in the various formalisms to follow. 
This model is based on the fact that the rotation group S0(3) i a nat- 


ural configuration space for the orientation of a rigid body. 


The next generalization of the formalism is to describe the 
nonrelativistic rigid body with the aid of the theory of connections on 
principal fibre bundles; the definitions of this theory are collected 
in an Appendix. At the same time, Chapter III makes the transition from 
the Lagrangian formalism to an inhomogeneous symplectic version. Although 
no new physical content is discerned (or expected!) this modern geometri- 
cal description of the (nonrelativistic) rigid body is useful to have 
(as a point of reference) for comparison purposes with the corresponding 
Newtonian and relativistic space-time descriptions. Also, development 
of the relativistic formalism can be guided by analogy with the non- 


relativistic version. 


Chapter IV is mainly a review of the homogeneous formulation 
of symplectic mechanics (time is included as a coordinate) and its 
application to the description of rigid test body motion in a curved 
space-time. Dynamics is now discussed on a higher level fee that (of 
the traditional Lagrangian (or Hamiltonian) in that a (pre) symplectic 
potential 1-form or presymplectic 2-form will be the basis for dynamics. 
These two approaches, not necessarily equivalent, have varying degrees 
of advantage depending on the situation to which they are applied, an 
example being as simple as the facility of choice of Lagrangian or sym- 


plectic form for a particular system. 
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The material of Chapter IV is used in Chapter V to obtain a 
symplectic formalism of constrained motion of a rigid body in space-time 
(the situation envisioned is the centre of mass of the body being 
constrained to a given worldline). While no new effect is expected in 
the Newtonian theory, we do expect something in relativity, based on the 
following analogy: It is known (84.6) that a relativistic spinning rigid 
test body freely moving in a curved space-time is subjected to a force 
which causes the body to deviate from geodesic motion (in a fashion 
analogous to that of a charged particle in a magnetic field). In anal- 
ogy to the situation in Newtonian mechanics (holding a top fixed at one 
point causes it to precess), we expect that, if this body is constrained 
to (say) a geodesic, this force would be manifested in a precession of 


the spin axis. 


Space limitations make it difficult to define all the mathema- 


tical terminology which is used throughout. Only key definitions are 
given in the text, while a major part of the notation is summarized in 


the Glossary at the end. 


The Einstein summation convention is assumed whenever an index 
is repeated in an upper and lower position within one term: 


d 
Ki = , xX , while the range of an index may be either explicit or 
fo} 


implicit; the range for the latter is to be taken from context. Greek 
indices usually refer to space (-time) components while Latin indices 
refer to general (triad; tetrad) frame components. Most often capital-— 
ized indices (A,B,I,A for example) assume the values 1 , 2 or 5 


while lower case indices (a,b,a,8, etc.) range over 0,1,2, g) 
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Capital Latin indices are always raised and lowered with the Kronecker 


delta: ce 


or 6 AR » and lower case Greek (Latin) indices are raised 
and lowered with a Lorentz (Minkowski) metric in the relativistic case 


only. 


The organization of the thesis is in the Chapter-Section for- 
mat. The numbering of equations is by section, so that 1.1(1) refers 
to equation (1) of 81.1; major definitions are also referred to by 
number. To facilitate cross-referencing, sections are further divided 
by lower case Latin letters; for example, l.la. The occasional footnote 


is indicated by a superscript number at the appropriate location. 
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CHAPTER IT 


Analysis of Rotational Motion of a Rigid Body 


§2.1 Introduction. 


Although the subject matter of this chapter is familiar to 
everyone in its traditional form (eg, see Goldstein [16]), the reformu- 
lation of rotational motion in terms of Lie groups may not be so famil- 
iar in spite of its inclusion into several of the recent books on 
mechanics ([42], and [52] for example). Since the end result of this 
chapter will be the model of rotational motion of a rigid body that we 
use throughout, we develop in some detail the modern differential-geo- 
metric formalism used by Arnold in his now classical work on mechanics 


on Lie groups [2] as well as subsequent treatments CRASS Ie ea 


The first section defines a rigid body and states some of the 
relevent classical results, leading up to Euler's equations. The differ- 
ential geometry that we use is standard; however, due to the great abun- 
dance of notations, §2.2 is devoted to presenting the notation, as well 
as some of the standard results (without proof) that will be used in 
the sequel. Useful references would be Kobayshi and Nomizu [24], Warner 
[59], Godbillon [15] and the surveys found in [19] and [54]. We discuss 
next the inhomogeneous formulation of canonical dynamics with emphasis 
on the Lagrangian viewpoint; this is then applied to a description GE the 


rotational motion of a rigid body. 
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§2.2 Review of the Traditional Formulation of the Nonrelativistic 


Rigid Body. 


(a) Foundations. Like that of a point particle, the concept of 
a rigid body is a mathematical idealization representing physical real- 
ity in an approximate sense. In the domain of Newtonian physics, this 
approximation is actually quite good so we assume the following defini- 


tion is valid: 


(1) Definition. A rigid body is a distribution of mass in space sub- 
ject to the constraints that the distances defined by each pair of 


points contained in the distribution shall be constant in time. 


The definition implies that at each moment of time the knowledge of 
four noncoplanar points of the body determines any other by the con- 
straints. '"Space'' referred to in (1) is Euclidean 3-dimensional space. 
From a philosophical viewpoint it may be desirable to begin with the 
more general "“affinely-rigid" body studied by Slawianowski [44]. This 
is a mass distribution in which the affine relations are constrained to 
be constant in time, no reference to metrical properties being made 
(physically, affine rigidity corresponds to the possibility of homogen- 
eous deformations of the body). However, we are mainly interested in 
rigid bodies as in (1). As we shall see, the definition of rigidity in 
relativity is much less straightforward and indeed, strictly speaking, 


can't be done at all. 


(b) Kinematical considerations. The remark after (1) implies 
that fixing a set of three linearly independent vectors at a point in a 


rigid body characterizes the configuration of that body. Thus, the 
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complete configuration of a rigid body in space, R? is specified by 


(i) choosing an orthonormal frame (K, sk, 5k) for. 


parallely translated to each point of Re ) - the space frame; 
(11) giving a point of the body (eg, centre of mass) ; 


(iii) giving the orientation, relative to the space frame, of an ortho- 


> > > 
normal frame (e, ey 0&3) fixed in the body at the point of (11). 


Formulated mathematically: xe Ry denotes the origin of the 
a 
body-fixed frame, (e, seen 2s) yand).the matrix (ey) represents 


the orientation and is defined by 
> 
(2) eneane, Kn 


It turns out that because of (i) and (iii), the matrix (e,) is orthogonal; 


that is, satisfies: 
(3) ay ha, Ge Ns 
This, together with (ii) and (iii) above implies that a rigid body has 


six degrees of freedom. 


In practical situations, one most often represents Ene sroOLa 
tional degrees of freedom by Euler's angles (see [16] for example) but 


this approach is of little help in our conceptual framework. 


Now, at any instant of time, t , we have seen that the 


orientation of a rigid body can be specified by an orthogonal matrix 
iY 


(a@junelteat one winstantse ta=et) the body frame was coincident with 
A fo) 
iG iE ik ’ ’ : 
the space frame Si emee (t) = one as time went on the orientation 
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: Ny 
changes continuously, to e, (t) » say. Thus, on a physical basis we do 
not allow space inversions of the body to occur; hence, the restriction 


that 
(4) cae eel 
A 


is reasonable. 


(5) Definition. The collection of all 3x3 matrices satisfying (3) 
and (4) forms a group called the special orthogonal group (in three 


dimensions) and is denoted by S0(3). 


Note that Euler's Theorem, that the most general displacement 
of a rigid body with one point fixed is a rotation about some axis, is 


just the statement that each a SO(3) has 1 as an eigenvalue. 


(c) Rotating coordinate systems. We assume that we are given a 
rigid body held fixed at one point and that its rotational configuration 


is described as in (b). 


Given a vector é in R? , it is useful to know how to relate 
> 
the expressions of the time rate of change of G in either the space 


frame or the body-fixed frame. To do this we first give the 


(6) Definition. Let t > e, (t) be a motion of the body and for each 


ate et (on(t)) denote the inverse to (ey (t)) . Then the vector 
I 
de 
(ocak Nhe) a tyes oe aD 
(7) WED E Y Oo, dette 


is the angular velocity of the body at time t 
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In terms of (7) then, 
(8) “—} = (— en ae 


(To show this, all that is necessary is to note that: 
7 Ne ‘ 
pady: Ey (G ee on kp = Cenace. kp and differentiate with 


respect to time using (7).) 


(d) Angular momentum and kinetic energy. As in (c), we assume 
that the rigid body is held fixed at one point. The physical quantities 
that we shall now introduce: kinetic energy, moment of inertia and 
angular momentum, are all defined with respect to this point ang, git 
general, depend on this point. These quantities also depend, in a funda- 
mental way, on the distribution of mass of the body. While a mass dis- 
tribution can be quite general as in the continuum theories of Noll and 
Truesdell (see [57] for example; a body is a measure space in this 
theory), the definitions can be carried through quite generally if we 
assume the body to be composed of a finite number of point mass parti- 


cles, {m, Opes warn NiMESay . 


As usual then, we define 
> > > 
(9) L = } i Cony 


as the angular momentum of the body (with respect to the fixed point); 


> + 

Ly and Vy are the position vector and the velocity vector respectively, 
of the a-th particle, Mm, Next, define the moment of inertia tensor 
as 


(10) i Gh AGG ae Alias Se ee 
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> 
where the XK are either the components of ty in the space frame, 
in which case I need not be a constant in time, or in the body frame 


where I is necessarily constant. The expression for angular momentum 


(9) can now be written as 
(11) (a= ee tbe Le = Tg, 


The total kinetic energy of all the mass points is 


(12) T = +) mv : 
or 

(13) T = Sued d 
or 

(14) Wes > Tet aie 


We call T the rotational kinetic energy. 


Suppose oe and = denote two points of the body. The 


— 
following formula relates the moment of inertia tensor I at x to 
1 


> 
the corresponding I at Xo» 3 
2 


2 
Tp = dap TDN Ogg & Nyy 
(15) 
ae > 
Yn = 4) = = 


(e) Euler's equations of motion for a rigid body held fixed at one 


point. The fundamental equation of motion for the total angular momentum 
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of a rigid body is 


Play 


(16) 


> 
where N is the total torque acting on the body, and all quantities 


> 
necessarily referred to an inertial (space) frame, (k,) . In our point- 


iy 
mass model of a rigid body, N is usually defined to be 


(17) ke ae 


> 
where Fy is the external force acting on the a-th particle (we do 


not include the forces of rigid constraint). 


Via equation (8), we can transform (16) to its expression in 


the body frame: 

& > > 

(18) a) = (ae ex L 
> 

dw) a ie 


19 I (— +w X Iw 
(19) Ge Mae 


i 
ZY 


(I is constant in the body frame!). These are Euler's equations. They 
take on a more familiar form if we write the components in the body frame 


ere dives’ 
Ce re x 


(20) Te eee eee ae Tee nye 


for example, if the principal axes of the moment of inertia tensor ([16], 


page 151) are chosen to be the body-fixed frame so that 


tL : sk: 7 a 
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I = diag (1, 1,514) 
then Euler's equations are 


“1 D3 gare 
iT Ww + (I1,-14)w w = N 


and the other two obtained by cyclically permuting the indices. 


§2.3 Notations of Differential Geometry and Lie Groups. 


(a) Manifolds. All manifolds will be assumed to be smooth (c ) 
paracompact, Hausdorff and connected. Such a manifold is often indica- 
ted by a capital Latin letter, M,N for example. The phrase ie 
local coordinates x’ " means that we have selected a chart (U cM, 
pate alae R') at some point p eM where M is m —- dimensional and x” 


(1 boene! < m) 41s the da-th coordinate function of x ; U > Re 


(b) Tangent bundle. If M is an m - dimensional manifold and 


x €M then TM denotes the tangent space to M He se ROU Sg TM 
x€M 


is then the tangent bundle of M with 71: IM>M, (xv) soe thee pro 
x x 
jection. TM and TM are the cotangent space at x , and cotangent 


bundle respectively. 


Let N bea n —- dimensional manifold and ¢6:M->N a 
smooth map. We denote the differential of @ , a vector bundle map, by 
k * x 
>, : IM + TIN and the dual of this map by $6 ! TN => ITM. In a local 


, at : ; 
coordinate system x” of M and y of N, the coordinate basis of 


* Ooo 
TM Lot 2 oo e = g yeeanaatoreat Mes {ax} socthat if v=v —~ Ee 
x a oO x O 


Ox Ox 






6 tol pepe oto lot ae 


} fel f7h7 8 
ce 
— 
7 Lsptap otagik v é OAD 
; nig 1B OA . 
is : ia = ¢€ = 
_ 
ID 2d wy 6) 
a2 g > ‘neds 
x 


, 





Ee Ag . press atiw M Jo ofhimed joegged Sis) 
ur . 4 
7 7 


S = 


‘\ 
late ie -coljgag 
i Py al 
be : 





-¥ Levis obee94 ov. Fee 
es a < 






_ 
= we eo 


i | FF diy + lis? is sor Lai 
> a 






1% »*  ygoa we hin: T doa By > , 16 





Pte gee, oe 


14. 


L 
A e 
TM then ,, (v) =v” a, where or are the components of @ and 
ox doy : 
8 % * 3o- a 
if = = —— 
fo} 0; dy € eee , then > (0) < Oo; dx 


(l)eeDeftinition., A map 80; Mo N is a submersion at x © M iff 


>| ; TM me TN has maximal rank of dim N < dim M 
x 


Let yeN. If at each x€E ay eS) » ¢ is a submersion, then on Gy) 


is an embedded submanifold of M ([59]). 


(2) Definition. A map ¢?: M-*+N is a diffeomorphism iff Nee > N>M 


exists and is smooth. 


CmebeciniticomeAitapmetsit = Newisean embedding iff 15. 3Mi7 1(M) cN 


is a diffeomorphism. 


(c) Vector fields and 1-forms. Recall that a vector field Eas 


a smooth section of TM ,&:M-+ TM such that 1°€ = id, and that a 
1) 


x 
1-form, o , is a section of TM. ~ Locally, E(x) = E*(x) and 


ox |x 


o(x) = 0, (x) dx” ¥ (often, we do not use “a explicitly). We denote 
x 
by %(M) the set of all vector fields on M and by % (M) , the set of 


all l-forms on M. 


(d) Integral curves. Let c : R+M be a smooth curve. Denote 


by 1: RR = RXR the inclusion t > (t,1) 


(4) Definition. The lift c' of c is amap c':R-*>™ defined as 


A ‘ QO 2 : 
Ce= tego Am Li local coordinates, if t > x (t)) is the’ curve in 
dx 


ac lhe ago) EN 
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Let c be a smooth curve in M and € a vector field on Me 


(5) Definition. c is an integral curve of € iff GlUn=""0Cc). Local-= 


a 
ly, this is just is = E(x (t)) : 


(e) Lie group. A Lie group G is a set with a group structure and 
a manifold structure compatible with the group operation (that is, the 
group operation is ee \eelr BGe @icha croup, n,¢,c Co then Chee )ae ne 
denotes the group operation; e or id will denote the identity element 


OL se ow. 


The Lie algebra g of a Lie group G is the tangent space to 


G at the identity: 
(6) gan 1G 


Let M bea manifold, G a Lie group. 


(7) Definition. A map ~: GxM>M is called a (left) action of G 
On mMe Let vila) 2M > MS x > w(a,x) is a diffeomorphism and 


p(ab) = p(a)ow(b), we) = idy. 


It is evident that G can act on itself in three natural 


ways. Let aeG be fixed and b any element of CG. 


(8) Definition. Left translation by a in G is the map de y Go Ges 
b > ab ; similarly, right translation is the map 1. ;+ G@ 7G , 
b > ba. An inner automorphism of G by a is the map 
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Note the following properties of the induced differential 


maps: 
Cee eae (ie a) eal Ge > TG 


so that 


(id) (he ak 


ee G4 : vie 
be Tob is an isomorphism 


The differential of the map i, “aa > .Ger, 


(i) x > T.6>+T _16 is quite important. Consider its operation on g 
aba 

clearly (Che : g>g . This evidently gives us a representation of G 

in Aut(g) , the set of all automorphisms of g . This representation 


is called the adjoint representation and denoted by 
(9) a Ad, = (i) : 


Let, & € x(G) be a vector field. 


(10) Definition. § is said to be left-invariant iif for each) ae, G., 


had = bok, eit Gye cee S Le - related to itself. 


If & is a left-invariant vector field then 
(11) @oR |) =e 
a’ * ie e 


go that this map {left-invariant vector fields} + TG =g is a Lie 
algebra isomorphism where the Lie bracket [&,n] = En-n& defines the 


1) 


Lie algebra structure on the left-invariant vector fields: 
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(12) Definition. o is a left-invariant 1-form (or Maurer-Cartan form) 


x 
TET rorea lle ea. oo Ga. (2) Cea Om, 


Letmxenme akon ibe aabasishote etawhereagn =jdimrGy 
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al n * * 
ikeeen Seo uedenotesstheldual basis for) g = 1.6 . There exist con- 
K 
stants CAR such that 
(13) bien) Sues 
A’oB AB Se 


and are called the structure constants of G with respect to the basis 


{X,} of g 


If A € g , denote the corresponding left-invariant vector 
field by A-e *(G) 3; similarly, 1f O€ 2" ,» denote the corresponding 
left-invariant form by Came Te then follows that TG is parallelizable 
by {X,4 IKEDA esd 8 similarly r'G S Goce » The Maurer-Cartan equa- 


tions are 


(14) dx = 5c onl a oe 


(f) Exponential map. 


(15) Definition. A homomorphism dé: R2+G is called a l-parameter 


subgroup Of Ga. 


Leto) exp, (A) be the unique l-parameter subgroup of G 


d 3 
such that Ti} =0 (exp, (A)) = A where A€Eg 


(16) Definition. The exponential map exp: g~G is defined by 
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exp(A) = exp, (1) . Further details and properties may be found 


dint) [59]. 


TethercesesthatmGmisea matrixecroup, exp is eiven by 


matrix exponentiation: if Aeg , then 


exp (A) =etAt+ aA + see 


where e€G again denotes the identity element. 


(g) Example: the general linear group. Clin)e= Chin, Reel seene 
group of all invertible linear transformations of Ree Rae ee 
(x? ; a=i1,...,n) be the standard coordinate system on Rk” . We iden- 
tify each ae GL(n) with ite matrix representative relative to the 
basis oe of! RX: let X= x8, eR ; then a(X) = ee a(d.) = 
> Caer} re According to this last expression, we can view a either 
as mapping the basis bg to the basis {a 9g) keeping the compon- 
ents with respect to each basis the same, or reciprocally, the basis 
We 


{a} remains fixed and the components eo change to a B X 


examine these actions in more detail in the Appendix. 


Left translation in GL(n) is accomplished by matrix multipli- 


cation: let ae GL(n) ; then & (bd) =a*b for any be GL(n) . In 


pe: Similarly, let v € 1, Gli(m) 3 then 


O = 
components, (2, (b)) B a b 


(24) = ave Tob ChiGn) 


The Lie algebra of CLG) , gv(n) = T,GL@) is diffeomorphic 


nxn : ; 5 
to R , the space of all nn matrices. The result of multiplica- 


tion in g%(n) of two elements A and 8B is their commutator: 
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[A,B] = AB-BA € g&(n) . We choose as standard basis for gi(n) the 


matrices (Ey) defined by 
(17) @;) = 6 


Denote by (eS) the standard coordinate functions on GL(n) and by 


(8) the inverse to (e,) Pee nen. 


(18) oe = en a. and 
J J oe, 
i 
(19) | a oO ae 
J O J 


§2.4 The Rotation Group in Three Dimensions. 

(a) In §2.2 it was argued that the group S0O(3) serves to 
describe the configuration manifold of a rigid body held fixed at one 
point. In this section we examine in more detail S0O(3) and its role 


in the description of a rigid body. 


Capital Greek and Latin indices both run and sum over 1, 2, 
and 3. We assume that a global space frame tk} is defined at each 
point of R? in accordance with Euclidean parallelism; {n°} denotes 


the dual basis. 


We recall from §2.1 that 
(1) EO CG eae ECL (CO) nas ee deta al) 


S0(3) is a closed Lie subgroup of GL(3) and we consider it as a 
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subset of GL(3) ; in particular, the usual coordinate functions ex > 


the coordinate vector fields eu and 1-forms de, will be used as for 
2) de, 


6L@3) a. 
(b) Angular velocity. Let 


(2) Ce Re sc) a, ee > e) (t) 


be a smooth curve in S0(3) , representing a motion of the body. The 


tangent vector to the curve, e(t) = ey, (t) a a= =) ,» can be 
de 
written as 
‘ Meet BA 
(3) e(t) = en or EB : 
showing that the quantity 
x Bala 
(r= Qy Ep 
(4) 
Been lee. D 
TN =e, oF 


is an element of the Lie algebra so(3) of S0O(3) and that 
(5) a(t) = 2. 


To see the physical meaning of % , recall that the body and 


space frames are related by 


(6) e, =e, & 
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which defines the orientation of the body at time t . Then 
(7) ene=t eee 4 


defining the (1,1) tensor 


=) 
il 
2 
(D 
8 
@D 


(8) 


Transforming to the space frame, we see that 


ie A DeN IP 
(9) Q = (8,2, en) ot IE | os 


defining the space components of @ : 
(10) ly = 6, 


Equation (7) is the statement that Q is the angular velocity 
tensor; (8) is its expression in the body frame while (9) is its expres-— 
sion in the space frame. The associated angular velocity vector 2.2 C7) 
is 

Q =e, =- s Celene (body frame) 


(11) 


(space frame) 


B= 

ll 
—) 
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il 
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To recapitulate: we have defined two tensors, Q from (3) 
and (4) and Q from (7) and (8) which lie in completely different 
spaces: € so(3) while Wis 7) @) . ® is the usual angular velocity 


tensor and has exactly the same components as { ; whence, in our model 
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of rotational motion (being S0O(3) ) is the angular velocity tensor 


in the body frame. 


(c) Lie algebra of S0O(3) . By differentiating the defining 


relations (1), one obtains 


(12) so(3) = {Ae g%(3) : Atay = 0} 


A suitable basis for so(3) is then {E,,=E)-E,,1< AcB< 3} 5 


more often, we shall employ the matrices 
1 
C13) E,=- mae E_, = -e E 
A , * 
Let iE} Mbe the dualibasis:of .g «% 


(d) We now give a more formal description of (c) in S0(3) 
Recalling the definitions in §2.3 and especially Example (g), equation 


(4) is 
(14) : oy Oe 


and if we identify © with eet nen 
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(a) Coordinates. Recall (§2.3e) that 1TS0(3) ~ S0(3) x s0(3) 
The left-invariant vector fields Mal area basis of T_S0(3) for 


each a € SO(3) , so that (e,,2°) forms a coordinate system for 


TSO(3) . We also have the coordinate expressions 
Ces 
16 i Mes 
a) oe as: cca 
5 
al NTE AK ls r 
E = -e on de 
(17) 1ff 
ee ak 
de, =-e € AL E 


(f) Maurer-Cartan equations. Relative to the basis {Et of 


so(3) the structure constants of S0O(3) are Cae = con so that the 


Maurer-Cartan equations are 


(18) dE = -F cf gBAg 


(g) Left-invariant metrics. A left-invariant positive definite 
metric I ona Lie group G is a (2,0) - tensor field such that I, 
is a positive definite symmetric bilinear form on TAG foneall ac. G 
and satisfies oy ie= 1... Note, that, this implies) that 1 is cbtain- 
able from a positive definite symmetric bilinear form Tig on g: 


* Pa (pe . 
L, = (2 a I;q- Denote Tia by I and consider I as an isomorphism 


x 
Pace so that 


kw * y 
(19) 1,0) = @ TR 1) 4(W)) eTG if ve T,¢ 
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The metric I determines a function L: TG > R defined by 


loo 


(20) L(a,v) = 7 1, (vv) : 


The following diagram is a useful summary: 


a a 

(21) i | | 5 
* 
a 


In our S0O(3) model of the rigid body, I is the moment of 
inertia tensor (necessarily with respect to the body frame), and has the 


form 


(22) Teg yee Hie al 


Now, & = of Ea is the angular velocity; the quantity 


~ A Baek 
(23) S = I(Q) = SA E = Tas 1B 
or (amounting to the same thing) 
2 ~A 
(24) Si) =)S, 8 
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is the spin angular momentum of the rigid body. This is justified by 


writing 2.2 (11) in terms of the body frame (e,) : L= Ly of = 


A 
B 
Taps? ae and noticing that the components La = Sy é 
The function defined by (19) is now written 
ai -i > 
(25) L(a,2) = 51,9 =F 1h, ta 


with the identification I+>+L and Q<+Q if there is no confusion. 
We recognize (25) as the usual expression 2.2 (14) for the rotational 


kinetic energy. 


§2.5 Canonical Lagrangian namics. 


(a) The modern geometrical formulation of saetiseg Lagrangian 
and Hamiltonian mechanics can be found in several books devoted to the 
subject (Abraham [1], Souriau [48]; Chernoff and Marsden ieee on; 
infinite-dimensional systems) as well as books on differential geometry 
({15], [51], for example). This section derives a Lagrangian formula- 
tion from the classical Euler-Lagrange equations, then sketches the 


general inhomogeneous formalism. 


(b) Let Q be the configuration space of some physical system 
with n degrees of freedom (so that dimQ=n ). TQ is then the 
state space (or velocity phase space) of the system. Locally, we have 

“ a Qa a 
@oordinates: (q.)»eon 0 —and 9(q iv), on) TQ. Let “iss fO°> & be 
a Lagrangian describing the dynamics of the system. Classically, a 
motion of the system is represented by a curve t * (x7 (t) sv (t)) in 


state space which satisfies 
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and the Euler-Lagrange equations 
d OL aL 
(2) rs Omere age O 
oq 


We first postulate that t > (x" (t),v"(t)) sean iuteeral 


curve of a vector field Xe #(TQ) . Then, 
a dq” a a 
(3) is = Redque=ax(q )) =. 
dig Ol Gio OL 

(4) at Stet = wa! ’ 
whence (2) is equivalent to 

Il; L 
(5) xR) dq” - =. dq” = 0 

ov oq 


Expressing everything in terms of Lie derivation by X , 
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L L 
dq”) = Ly A) dg® + — Ly(dq’) and 
ov OV 
L, (dq”) = d(Xt#dq’) = dv” 
and defining the 1-form 


(6) UF eo Je 83 


equation (5) becomes: 
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We define two more quantities; 
(8) Wo do, ‘ 
(9) E = abe, = - L 
O 
OV 
Since L, (6, ) = d(%6, ) +X do, = Sart v°) - xa Ww, then we obtain 
the final form: 
(10) X 3 a, dE = 0 
9-1. 
If L is non-degenerate (that is saa is nonsingular on TQ) then 
dV OV 


w. is a sympletic Porm 2 


L E is the energy of the system. 


The steps leading to (10) are reversible, showing that the 
classical formulation (1), (2) is equivalent to the symplectic formula- 


tion (10). 


(c) Canonical mechanics. Again Q denotes configuration space, 
* 
TQ state space, TQ (momentum) phase space and L: 1Q > Rema segtan= 


x 
gian. Let (a sp,) be a local coordinate system on T Q 


(11) Theorem. ([51] pg. 143). There is a 1-form 05 € x (7 Q) defined 
by the property: for each u € TQ and for each X € Tigo we have 

<X, 60 ier <T1,X,u> where T : TQ +Q is the bundle projection and <,> 
is the natural pairing between ron and et in one case and TQ 


* 
and TQ in the second. The 2-form Wo = - dé. has maximal rank and 
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Bs 
thus defines a symplectic structure on TQ. In coordinates (q",p,) 


on r'Q OF “4p dq* and w = dq” A dp 
we eG a. Oo : 


Let Lec (TQ) and L. denote L 
q TQ 
q 
x 
(12) Definition. The fibre derivative of L is amap FL: TQ*>TQ 


defined by: for each u = (q,v) € TQ , FL(u) = (q,p) where 


* 
pe Bag has the following expression: 


p(w) = Sk ce forall vw & ne 
In a local coordinate system (q°) > v= v ~.. and FL(v) = oe aq F 
oq Ov 


* 
(13) Definition. L is a regular Lagrangian iff FL: TQ*>TQ isa 


regular map on 1Q (ie, dFL is nowhere zero on TQ ). 


(14) Proposition ({[1], pg. 116). L is a regular Lagrangian iff FL 


x 
is a local diffeomorphism iff FL (ow) is a symplectic form on TQ. 


We remark that, in general, the fibre derivative FL is not 
linear as a map on the fibres of TQ. If ve Hay is the velocity, 
then p = FL(v) is the canonical (or conjugate) momentum associated to 
v. The symplectic manifold (T°Q,W,) is the arena of Hamiltonian 
mechanics while (TQ, 0, = FL (w.)) is that of canonical Lagrangian 


dynamics. 
We need two further definitions: 


(15) Definitions. Let L be a regular Lagrangian. The function 


AM LOur Re RAG. Vv) “aay, bL(v)= is the action of the Lagrangian 
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L. The function E = A-L is the energy Ofte ln, 8 lne coordinates 


(q” TD on TQ , 


(16) a= oh yO 
Ov 
Ov 


Let L be a regular Lagrangian. The symplectic form 
w= FL (w,) on TQ has maximal rank and thus provides a linear iso- 
morphism between T 1Q and TQ for eachasuse TOs, Denote by Ww 
the induced isomorphism between #(TQ) and ¥ (1) SOM ciated ke ees 
any vector field on TQ then w(X) is a 1l-form on Q. The canonical 


Euler-Lagrange equations are 

(18) X = 0 (dE) 
or equivalently, 

(19) X 3 w= dE 


which is seen to be the same as in (10). The motions of the system are 
represented by integral curves in TQ of the vector field xX» satisiy— 
ing (19). It is shown in ([1], pe. 120) that X defined by (8) 7isea 


second-order equation on Q ; 


(20) Definition. A second-order equation on a manifold M is a vector 
Eieldse\euoheeiMe SucHythat Tos = idiy where 7 ; IM > M and 


Tee ee eI Ls its differential. In a local coordinate system 
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(21) xX = a a + £"(x,v) ae 


Last, we remark that the Lagrangian structure on TQ induced 
by the canonical structure on TQ is not the only way a Lagrangian 
formalism can be formulated. An alternative approach, making use of 
the so-called almost tangent structure on TQ , is used by [23] and [58] 


and discussed briefly in 85.3. 


§2.6 Canonical Dynamics of Rotational Motion of a Rigid Body. 


(a) The canonical Lagrangian formalism described in the last sec- 
tion is now applied to the configuration space S0(3) representing 


rotational motion of a rigid body as formulated in 82.4. Henceforth: 


Q = S0(3) 
TQ = SO(3) x so(3) 


De elOr> Ree (d,s) = $1 (25)) = 51 2 2 


q 
(1) 


I 
n 
i 
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eo) 


FL : TS0(3) + T S0(3) , FL(Q), 


It 
ip) 


9 *g" Q 
eek (T'S0G)) 4 = 3, 


* 
o. = FL (6) = Tp E 


The canonical symplectic form w, = -d6, can be computed using the 
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Maurer-Cartan equations 2.4(18): 
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(2) w= Lap HCehEL ee 5 Sa_ EAE 
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where we still write S, =I ee for convenience and 


(3) S =e S 


The energy (see 2.5 (17)) is 


if B 
(4) E= 35 lap ne 
so that 
(5) aig Sry ee : 


Let X € ¥#(TSO(3)) be any vector field, say, 


(6) x= PAp + A) 
Phy) 
Then 
B B. cA By 2A 
Cia Xa wy = ~(S,pP + Tape JE + Tape di? 


so that equation 2.5 (19), X-+ w = dE , implies 


(8) pv = and 


py 


eee ) 


) is the inverse to ap 


Note that X given by (6), (8) and (9) is of the form 2.5 (20) ,0 a 


second-order equation. 


eul 
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(c) Integral curves. Let t > (e, (t), A(t) be an integral curve 


of X ; physically, this integral curve represents the evolution of the 


system with respect to time t . Then, 
fy 
de 
A is i Re Sb, Ceca 
(10) AG X de, XA( e, € ALE ) EAR Q2 eo 


where we have used 2.4 (17). Equation (10) is merely a regurgitation 


of the definition 2.4 (7). 
Last, for angular velocity we have 


ash A ~AB C 
destabedion, wal Ee = 
i X dQ aL Sac QQ 


(11) 


which we recognize as Euler's equations, 2.2 (20). 
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CHAPTER ILT 


A Canonical Lagrangian Formalism for a Nonrelativistic Rigid Body 


§3.1 Introduction. 


We saw in Chapter I how the motion of a rigid body held fixed 
at one point can be modeled as a canonical Lagrangian system on the Lie 
group SO(3) . The generalization to include motion in space according 
to a fairly general Lagrangian is made in this Chapter, partly for its 
intrinsic interest, partly to develop more formalism as preparation for 


the space-time formalisms which will be met in Chapters IV and V. 


In simple situations (for example, a gyroscope falling in a 
gravitational field) the rotational and translational motions occur 
independently and one takes R> x $0(3) as configuration space with a 
Lagrangian such as: L : T(R?xS0(3)) > Roe meee su) ace > mv + Lee Ge ~ 
mo(x) , m being the mass of the body. However, we shall generalize in 
two directions: First, space shall be an arbitrary 3-dimensional Rieman- 
nian manifold = , representing the configuration manifold of a point of 
reference fixed in the body (usually the centre of mass). The physical 
significance of allowing *% to be so general is not clear; however, it 
may have some interpretation in the mechanics of generalized continua. 
Our main purpose though, in the generality of 2% , is having the present 
formalism as a point for comparison with the later space-time theories. 
The second generalization has already been mentioned: a Lagrangian which 


can include electromagnetic-type interactions and spin-orbit like inter- 


actions as special cases. 
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§3.2 Bundle of Orthonormal Triads as Configuration Space. 


(a) It was remarked earlier (82.1) how the specification of a 
body fixed triad of vectors (with respect to a set of space axes having 
the same origin - more about this later) leads to the description of 
the rotational motion of a rigid body by S0O(3) . Thus our representa- 


tion of the general configuration of the body consists of specifying 


(i) a point x € y , a 3-dimensional Riemannian manifold with 


metric tensor y : (centre of mass). 


(ii) an orthonormal frame e = (e,) of T,, y : (orientation with 


respect to centre of mass). 


The totality of all orthonormal frames at all points in 2» 
leads to the choice of SO(2) , a principal fibre bundle over 2 asian 
SO(3) as structure group, as configuration space for the rigid body 
(see the Appendix for the definitions, or [24] for a complete discussion 


of principal bundles). 


(b) SO(Z) : definitions. SO(z) is a sub-bundle of L(y) and 


is characterized by the metric on 4: if (x,e) € L (2) , then 


(1) (x Ve) ~e SOCx) 1th Y,-6€4 2p? = 5 AR 


It is clear that if ae $0(3) then 


(2) e € SOF iff eae SO 2% : 
x x 


I 2 


and this shows that S02 S0(3) 
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This process of obtaining SO(z) from Les called 
reducing L(Z) to SO(Z) ; SO(z) is one of several examples of G - 


structures that we will be meeting (cf: [51] for the general theory). 


A local coordinate system on a neighbourhood U of 
uot xse)mces0(.)e willeatili abe denoted |by (x! e,) (see the Appendix) 


but with e = (e,) subject to the results of §2.4. 


The linear connection JI. on © that we will use is the 


ordinary Levi-Civita connection of the metric y: 


ie Hh oye ee 


0 
AA 


The so(3) - valued connection 1-form is 


W = sw EAR 
ae eee | We i a FS 
(4) a @ (de, + Le, e, dx ) 
‘Gaess 
E aes Eee Saat 


The canonical 1-form is essentially the same as AG is 


Se (Ge) 


(5) 


at a point (x,e) € SO@) 


We take fe, EI as a basis for ¥(SO(2%)) with 
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and the corresponding dual basis is denoted by (on Eo 


(c) The interpretation of the various quantities is essentially the 
same as in §2.4. on are the components of the angular velocity in the 
body frame e at a point (x,e) € SO(Z) when the centre of mass is at 
x while ee ea projects (under 7, : TSO )=2sIL 9) cto the velocity of 


the centre of mass: Ty (ve) = vie —F Se} 


Ox % 
Orientation of the body has been defined with respect to the 


coordinate basis i} (at least locally). In a general setting it is 
Ox 


difficult to determine a priori whether the body is rotating as the 

centre of mass moves through space since I #0. After the path c in 
E of the centre of mass has been determined, one can compare the behav- 
ior of the motion t > e(t) to a chosen frame field which is parallely 


propagated along c , at least in the present situation. 


(d) We list for reference the structure equations of the connec- 


tion (see A(15) and A(16)) in this case: 


(8) do = -<,, EA 8 
jae eee 07 oC Seen EG. ek all 
(9) dE = = mE pce AK -F 5 R rai AQ) 
where sR meee = oe is the so(3) - valued curvature 2-form of 


the connection. 
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§3.3 Canonical Dynamics of the General Motion. 


(a) The programme of canonical Lagrangian mechanics outlined in 
§2.5 is now applied to the configuration space Q = SO(%) . The Rieman- 


2 2 e AN 
nian metric on % is written as either y = lee dx ® ae or 


A B 


Y= OAR @ ® 6 . We assume the rotational motion is described as in 
B25, 0.. 
The Lagrangian) = > TS0(2) > Rk is 
cee A B il A.B AUB 
(1) Le 7m 6 apY vo + 5 Tap” QO + Kav Q 
A 
+ A,V + Bt - md 
where TAR is the moment of inertia tensor with respect to the centre 
of mass, 
(2) T = a 6 ene : translational kinetic energy 
trans 2 AB 
(3) cc Se Ag® : rotational kinetic energy 
rot 2 AB 


and (2) and (3) are both defined with respect to the centre of mass. The 


term Kyav » with Kap depending on x only, represents a sort of 
"spin-orbit" coupling while Aw and Bt are velocity-dependent 
potentials (A, eraticl By are functions of x only). Finally, m¢ 


is the potential energy and again ® dis a function of | x only. 


(b) Special cases. Unfortunately, the Lagrangian Gi) es) too gener= 
al to generate the corresponding equations of motion in a reasonable 


amount of space. To illustrate the method and typical results, we exam- 
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ine two special cases; 

Case (i): 
(4) L = =m 5 apY V ace ar + K ve 
where K is constant, and is "small" compared to nag A 
Case (ii): 


ol PN 1B ik B A A 
(5) L = zm O apy vo o+ 5 Type 2 +m AV +m BAS - md 


where A,B and @ are functions depending on x € 2 only. 


(c) Canonical 1-form for the Lagrangian (1). Define the quanti- 


ties 
(6) A ie Hie kee aa 
Pe ag atau | eet 
2. B B ¢ 
(7) Vee TAN e PK vent Bea: 


these are the momenta conjugate to v, and aN respectively. Then, 


A ~A 
(8) a TOs ar LAE 
is the canonical 1-form of the Lagrangian (1). 


(d) ‘Case (i). We first compute the energy & associated to the 


Lagrangian L given by (4): 


(9) E = vO, ais eae 
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so that 
(10) dE = (mv,tk, 0 dv" + (Ly, +K, v) a0 
From the 1-form (8) we deduce the symplectic form w = -d0_ , 
(11) w= -m dv,r0™ + K,,0'AdQ + €., 0% EPAe” 

= Ip ange + Kea fleas 

Pee Oe pec ao 

Let 
(12) Z = vee, + PPE, +E" + oA 2 
ov 082 


be a typicalevector fieldton = %1S0(~)) © Using’ (1) and (12)> "the 1-form 


Z4W is found to be 


B iL 
(13) ZAW = ~ (Kp? : mE, ‘ EaBce me > “tune AB 


B BC. ~A 
~(Iyp + Kya” ~ €gpol oie ee 


A B A B B 
+(mv + Kap jdv + (KY + Tape Wales 


Then the canonical form of the Euler-Lagrange equations (2.5 


(19)) : Zi3w - dE = O implies 
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B B BC ree, 
(15) a Ky Ray i ce See 50 
(16) mv, + K of =nvy, + K ae 
A AB A AB 
B B_ B B 
(17) Keay + Tape = Tape + Noa V 


using (10) and (13). 


Equations (16) and (17) may be rewritten as 
(18) m(v,-v,) + K (oer = 0 and 
A A AB 
(19) Kya@r-v?) + Typ) = 0; 
then the substitution of (18) into (19) gives 


aie = 
(20) copes, ) (Paxti sh ORs 








Tf | 1K] | < vm||I|/| in any suitable matrix norm | ‘ | seLnen 
mi - KK has an inverse, which we denote by I . 
We conclude from (20) that 
Can) oe = oe and (18) gives 
(22) Me = ee 
To save some writing, set 
(23) Ry = eye - = ery ener: and 
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, BAG BC 
€ ABC QR? + € Ov 


cous ABC 


Pe 
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Substituting (14) into (15) and using (21) to (24) gives 


A _ =AB, = G 
(25) o = T'(@R, - R°K,,) and 


A -1,A —BC— ie D 
2 = = 
(26) a m ~(R -mK ph Ra + Kapt Topk ) 


where again a denotes the inverse to 


We have essentially solved the problem, the components of the 
vector field Z being given by (21), (22), (25) and (26). What is of 
interest, however, is the motion of the centre of mass in space, the 
time development of angular velocity, etc. This is achieved by projec-— 
tine Z into SO(Z) and 5 . Todo this let t denote the parameter 


of the integral curve of Z 3; physically, t represents time. 


The tangent vector of the trajectory of centre of mass is 


IR 


ee Che Ura h (el ne 
C27) Ste era ZAdx Z e,0 ve, 
From (22) we see that 
(28) ust = vs e 4 
Similarly, 
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= Z1(-e ean te - ve ve *0 iy Gap eete. 
(29) oe ce ate 
where ''e" = v'Vp is covariant differentiation along the trajectory 
Caer Ce) 
Next, 
(30) aut” = 2449 = 08 


so that computing of with the aid of (23) and (24) we get 
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For the linear acceleration, we have 
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or using (28); (29) and (26) 


T aN. ete 
(32) v= =e, € BC’ Q 
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(e) Subcases of Case (i) include: free motion of the body where 


Kap is peten to be zero in (d); (31) and (32) become respectively 
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(33) agree —]. SRC Q2 
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(34) vo =m 5 Shs Vv 


The other subcase is where 
(35) K = ime ie r 


a being a constant vector (in the body frame) which labels a point of 


the body other than the centre of mass. It then turns out that 


=-1 A iL 
m 


lee : , : : 
TIT = I-m KK is the moment of inertia tensor with respect to the 


point eh ee + i 


of the rigid body (compare formula 2 20l3)) 2 Lhus 
(d) with (35) is the description of the motion of the rigid body with 


respect to a point of the body labelled by n relative to the centre 


of mass. 
ee ee AB, 1 ACB A A. 
(f) CasenGi): L = 7m op” vt Ih ea Q +m A.V +m BAS mo 
is the Lagrangian (5). In succession we compute that 
(36) On = mv, + mA, 
(37) y, = S, + mB, 
Ih A ah B 
(38) B=fmv vit 5 Tye + mp 
A A 

dE =m v dv ae Sy age + mp 0 

(39) 
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so that 
W = -d6, seOrg using (30) to (41) 
A = Bin G A bon D 
= A - A 
(42) () = mv € ARCH fs) mdv , Q- + mle] 8° AG 
1 A eB oC 1 ea bela ie AB Kalle 
Sora Eapce *E + > S,p(E AE +> R Tae 
A.~B A,~B 
cs A 
aN di? AE OB © E 
If again Z= wee ote pg ar eee ave denotes a typical vector 
~A aH A A 
Ov ee. 
field in ¥*¥(TSO(S)) , 
Bac B B 1 MN Ko nA 
= — - = Vieie 
(43) Z4.W) (mé ,tm €apce vo+ Booey + Eee 5 S RunKA ) 
B B B.C B Bo Geek 
~ - - - 9) 
(In? + S ape M € apc Vv mV Bala M € RG B )E 
A B A 
= mV, ahi se Tape d%2 
The equation ZW = dE holds. iL 
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As an example, we can write the equations determining an integral curve 


of Z in three-dimensional vector notation for the special case that 


i \ aa ; > > 
Y= Sn, axBee “dx and B is constant: with B = V ~ A: 
; + 
(48) v = vxB - Vo 
(49) 8 = ikS + OB 
P ie > : ; : 
where v= (v') and S = (S,) . It is evident that (48) is an analogue 


of the usual Lorentz-force, and that (49) represents some sort of inter- 


> 
action of spin with an external magnetic field B. 
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CHAPTER IV 


Motion in Curved Space-Time 


84.1 Introduction. 


So far, our discussions of the rigid body have taken the 

form of describing the evolution of the relevant physical quantities 
such as position of centre of mass and rotational configuration, depen- 
dent upon a (fixed) parameter, the time. As the title of this chapter 
indicates, we now generalize the formalism to include time as one of the 
coordinates; the motion will then appear as a (unparameterized) curve in 
a suitable (extended) configuration space. In doing so, we can not only 
discuss the relativistic rigid test body, but also include naturally a 


parallel description based on Galilean relativity. 


Although we still rely on a definition of rigidity, such as 
2.2 (1), in a nonrelativistic setting, the corresponding definition in 
general relativity is not immediately obvious (an extended body cannot 
be considered at one time, since simultaneity can't be defined over a 
finite region of space). Any "reasonable" rigidity condition, such as 
Born rigidity, for an extended body seems to be too restricting ([41]). 
However, we may retain some of the properties of a Newtonian rigid body 
if we adopt Dixon's dynamical rigidity condition ([13]) which essential- 
ly Se that the classical relation (for: example, 2.4 (23)) between spin 
angular momentum and angular velocity is valid for the body. In this 
case, there exists a mass constant, m , such that the rest energy, Mig 


of the body has the expression 
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Ch) Mas itetol ea: + ( potential energy of interaction ) 
with applied fields 


We place an additional assumption on the dynamically rigid 
body, namely that quadrupole and higher multipole interactions with 
the gravitational field are negligible. Then the last term in (1) does 


not appear when there are no electromagnetic fields present. 


In the relativistic case then, we confine our discussion to a 
rigid test body which is assumed to be characterized by the quantities 
defined in §2.4. Our model (in both the relativistic and Newtonian cases) 
consists of an orthonormal triad of space-like vectors fixed in the body, 
determining the orientation of the whole body at a given point in space- 
time. Hence we think of a rigid body existing, at a given event, only 
in the tangent space to space-time at this event. In practice, the 


model represents, for example, a gyroscope orbiting the earth. 


With these understandings, we begin in the first three sections 
by reviewing the Galilei and Lorentz groups and algebras, Galilei and 
Lorentz structures on space-time, and the general programme of symplectic 
mechanics, which will be the basis of the model. Here, we are mainly 
interested in establishing the notation. Discussions of the Lorentz 
group can be found almost everywhere while the Galilei group and Galilean 
relativity is extensively reviewed in Levy-Leblond [30] and treated in 
several recent books [42], [52]. Galilei structures and corresponding 
comparisons of Galilean and Lorentz physics can be found in Havas [18] 


(traditional), Trautman [55] and Kunzle [28], for example. 
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The last two sections motivate and describe the evolution 
space, presymplectic form, its kernel, and the Beatie description of 
the motion. The treatment is, for the most part, that in Kunzle [28], 
although the Galilean treatment is new. We refrain from discussing sym- 
metries and conserved quantities since this is a deep topic in itself, 
and would lead us too far astray to give an adequate treatment here. 


(However, see [26] and [28], for example.) 


§4.2 Galilei and Lorentz Groups and Their Algebras. 


(a) Galilei group. Let (t,x) denote a point of Ro with 


oe pate a A we Sg 
x = (x ,x ,x ) = (x) . The proper Galilei group, Cale. scanepe 


defined as the group of affine endomorphisms of Ro of the form 
(t,x) > (t',x') = (ttb, Rxt+vtta) . 
(1) 


je) Se ee es R? y Ren oOC3) 


If v=0, we have essentially the Euclidean group of R? acting on the 
(affine) subspace t = const. For b=0, a= Q , and R= iid {7 we 
obtain the pure Galilean tranformations which are interpreted as switch- 
ing on a motion with velocity v , effected at time t = 0 in the above 
parameterization. Corresponding to a = OF and 9b8=.0)  isethe  bomogen- 


eous Galilerootoup,, Gale. 


Instead of considering Gal as an endomorphism group of Rt ‘ 


we will mainly think of Gal as operating on Galilean frames (to be 


made precise shortly) where (t,x) , (t',x') are the coordinates of 
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the same point of Ra in two Galilean frames. 
(b) Lorentz group. Here we consider the connected component of 
the identity only, and denote it by Lor . We can represent an element 


Eves Lore oy the matrix 


(2) 


Ae S0(3) AE, eR ; 


and where exp denotes the exponential map to a Lie group from its 


algebra (cf. 82.3). 


(c) For comparison purposes, we represent Gal and Lor in 


GL(4)) acting on the left on Rt : 


ae Gal<> a 


N 
a 
ov 


0 
R 
Le Lor< L= a Soe 


(3) b R+ CS) va’s) 
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(d) Lie algebras. Throughout the sequel we define 5) 

a 
( 0 ; Galilei case 

(4) ve ¢ 
pele, Lorentz case 
N 

: 
Let {EL} denote the standard basis of g1(4) . Take, then, as bases 


of the Lorentz algebra, lor, and Galilei algebra, gal, the matrices 


ae 
a IN B 
(5) i ees ay Le es BS 8 
Oo 
Pia eae 


Noting that LE generates so(3) , hence 


Peete 
{Ey = 5 an 


as a basis of gal (v=0) or lor (v=1) 


: eed! also generates so(3) , we often will use ee 


§4.3 Galilei and Lorentz Structures on Space-Time. 


(a) Action of Gal on frames of a 4-dimensional vector space 


Vie Let 
( (e.) be a frame (basis) of V , and 
(1) | (97) ‘the dual frame of v" SVin(V, 8), So. Enat 
en(e,) =o - 
Lét x 6.V ., x= 2 . Then, if Ge Gal , the coordinate transformation 
x? + x? = Ge au corresponds to a basis transformation arr en che of 


oe | i. ie * 







26 pied? .oar1 tj) lo: Yoveie@nd PMalyumge oils stom ‘ a 
tutes ol? ,Tse .padegle SebiTeqhst ,:0) i aaa Saye 


ry 
ty 
ae 
sy 
’ 
oe 
% 
Ms = 


iis ] if 2938 IOGE i a gade petzon 
liw wsgio. s (Oe Het ebogrs4 (ugha: Nag =: i “3 


b=) wet” me ow} Teg? Sa vtad: 





sei T-s9 698 "me (30 tea 
T8958 eae 2 Rae Lee al 


e3A0 miscy Jenéleacsa!'b-) q ip 2ome v4 aly iss) ie. SSizsn - oe 


hy py MO, co hae by aD: 
deild eas 1 
‘7 ne 





V , and the transformation of the dual basis of Wwe 6? = Gb a> 


Recalling 4.2 (3) we have 


G i = 1] G N =) 
eA Aguwers 

(2) cue b Cinta be, 
WHER yy IO re Slee 


so that the basis transformations are written as: 


Sens SN 
aa) Te e,@ ) . b 
> o 23.8 
ay CS Te 

(3) 5° = 6° 
BY Se R's 9 


These relations (3) imply that the tensors 


LAS ea 
p= 90 hela! 90 6 en ® ep 


are invariant under the action of Gal ; conversely, given a 1-forn, 
x p , 
y eV, and a symmetric rank-3 contravariant tensor, Y € V @V =, such 


that yy) = 0 , we can characterize an action of Gal on V 


In the following remarks, A = Ge) denotes a generic ele- 


ment of either Gal or Lor. 


(b) Bundle of {Gal | Lor} frames. We have defined the bundle 
of linear frames in the Appendix. The strategem now is to restrict our= 


selves to those frames (at a point) related by a {Gal | Lor} transfor- 
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mation (compare 83.26); this we make precise. 


Leta.y. be a ee four dimensional manifold (space-time). 


Suppose that at each point x € V we are given a frame e = 2 LV 2 
Oo 


e 
ao 
and that this assignment, x > e (x) , is smooth. Let G denote either 


of the groups Gal or Lor 


Forseach) x e€ Vo, define cal by: 
(4) Cam (ee) soaCe Vermette Cera ae Ne for some AEG. 
a x geo a 


Then set G(V) = vu GV e G(y) isea principalyfibre bundle over space- 
xeV 


time V with group G = {Gal | Lor} . Note that G(V) depends on the 
choice of x > e (x) . V together with such a reduction of L(V) to 


G(V) is called a G manifold. 


(c) Tensor characterizations of G(V) . Now we must separate the 
two cases according to whether G = Gal or G = Lor . Consider Lon) 


as in (b) and denote by Nab Ene gmatrix, ediag( =i) cn Deri neg: lic 


B 


smooth symmetric tensor field on V g = B48 dx” @ dx by 


a b 
(5) Gere eects) 


BoB 


a * 
where, as usual, g° = ep dx” denotes the frame of wy dual to e 


6 a b 


Clearly, this tensor field, g = Exe dx” @ dx” = n v) is invar- 
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Recalling the results of (a), we can similarly define a tensor 
characterization of Gal(V) . Assume that we are given a Galilean struc- 
ture on V as in (b). Define the following l-form, y = vy dx , and 


symmetric rank-3 tensor field, y = phe oy ®@ og 


ee 
oa 


fy, = 8° Go) 
(6) d 
| Bad = shoo 8 g560 


Conversely, suppose that we are given a l-form j = by dx” anda sym- 


metric rank-3 tensor field y = yor 3 @ 99 such that yor Ve = 0; 
define a Galilean structure on V by: 
a O 
e by = 6. and 
(x,e,) € Gal(V) <=> 
BY © leh ON) OG Sel Nee) 
oF Y Oe = 6.0 oe 


(d) Terminology. Let V , G(V) be as above, and suppose xeV. 


Q 


ASVECLOL © oe= ac or € ee is called 


vy, & #0, Gal 
timelike iff 
a-B 
Bug ems fe 8) iver 
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by eo =n0 , Gal 


spacelike pliant 
a-B 
Byg & & PaO a uLOk 


and, in the Lorentz case only, 
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Note that in the Galilei case  W defines a 3-dimensional 
distribution, x -> S.. c TLV , defined by: 6&6 € S ie oc Nii) S200 
thateis; S.. consists of all the spacelike vectors of TV eget exec =e 
iseintecrapies(1.¢,.. the lie bracket, [€,4] 5; of any two vectors, € , 
Qe ke S.. is again in > ), then each of the submanifolds {x ; x =t} 
carries an induced Riemannian metric (,) : if € = Ee e, and 


ia 6. © S. » e, being any Gal frame at x , then (€&,n) = 


We will be concerned with only those Galilean manifolds, 


(V,v,Y) , which are first order Piatoce that is, those for which Ww 


satisfies 


(7) dy = 0 


(e) Lor-connections. We do not need to say much about the connec- 
tion which we use on Lor(V) : it is merely the Levi-Civita connection 


induced by the metric g . The lor-valued connection 1-form is 


A A 


EA ane de Fy » where 
AC 
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w= E 


e el dx” 


ene ae ) 


Peel a 
E =- 7 € s) (de, + T 
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B 
fo) 
Se A,, oO eye vores 

6 (de. + Pay en (apie) 


er the By are the usual Christoffel symbols. 


(f) Gal-connections. The situation now is a bit more arbitrary 
than that in the Lorentz case. (y,W) does not determine a symmetric 


torsion-free Galilei connection uniquely, but it turns out that Cj28)) 
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feet) Weyyeeisea rirse order flat (f.0.£.) Galilei manifold, any con- 
nection is torsion-free and the set of all connections on V is ina 


1-1 correspondence with the set of all 2-forms on V. 


Ined £20.£. Galileilimanifold, (V.W,y) , we have, dj = 0 , so 
thata(locally)eave=—t 2aVereRa? shencé we can, find (locally) an adapted 


coordinate system such that 


( wae 
vy - oF 
| = ut a 
Bored @ 
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The components, in this coordinate system, of the most general Galilei 


connection are 
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If, in addition, we require that 


(12) vo toRY] oO, 


(85)p © 
then there are enough algebraic restrictions on R avg to reduce the 
number of independent components to 20 as for the Lorentz case. Equa- 
tlonecl2) .inesract, is -equtvalent to saying that onva flo.f: Calilei 


manifold the connection satisfies 
(13) dk = 0 


in the notation of (9) and (10). (This is expressed more generally in 
[28] than we need here.) Such a connection has been called a Newtonian 


connection ([28]). 


(g) General formulae for reference. Here V is any Lorentz mani- 


foldg(V,@eeOtet.O.t..caliler manifold (V,V,Y) with Newtonian connec- 
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rion eAsmusual.« ve= OM refers to the Galileigcase, Vv = 1. to the 


Lorentz, while G refers to one of Gal or Lor , g to gal or lor 
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84.4 The General Programme of Symplectic Mechanics. 


(a) In view of the increasing popularity of symplectic mechanics 
and the number of diverse forms the theory is being developed in, it may 
well be useful to include a brief survey of the field and to make explicit 
the form that we use in the sequel. Without going into detail, we hope to 
resolve some of the ambiguity of the varying terminology and usage of the 


various structures found in the literature. 


(b) Inhomogeneous canonical dynamics. This formalism (developed in 
[15], [51], and especially Abraham [1]) was discussed in some detail in 
§2.5. Suppose Q (a C n - dimensional manifold) is the configuration 
space of some system; we saw that TQ has a canonical symplectic struc- 
ture, Wo = sich » where oe is the canonical: form on TQ . M= TQ 
is known as the momentum phase space or simply phase space. A vector 
field on M Xe X#(M) is called locally Hamiltonian iff Ly w= 0 
(X generates a symmetry of Wo .) Then locally there is a function H 


such that X = x, where X. is defined by ee De dHee Hi eis*cealled 


H 
a (local) Hamiltonian for x, . Hamilton's equations of motion are essen- 


tially the equations for the integral curves of xy : 


In §2.5 we were given a suitable function L: TQ>R, the 
Lagrangian of a mechanical system. TQ is the velocity phase space or 
just state space. In terms of L , we defined a (regular) map 
Pines LQ TQ (the fibre derivative) and an energy function E : TQ > Res 
Via FL , the symplectic structure is transferred to 1Q , and the canon- 
ical form of Lagrange's equations obtained as integral curves of the 


vector field defined by xX, w = dE . The projections of these curves 
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to Q coincide with those obtained from Hamilton's equations when 
x a 
He: TsO)" R “is derived from E “via” FL “(that is, 4H = Eo (FL) ts at 


FL is a diffeomorphism. 


(c) Contact Manifolds. An odd-dimensional manifold M_ on which 
is defined a closed 2-form w of maximal rank is called a contact mani- 
fold and w a contact merverurerse Contac’ structures are basic to time- 
dependent mechanics and can describe either Hamiltonian systems or Lagran- 
gian systems. We give a simple example of the former - Lagrangian systems 


are similar. 


Let (M,w) be a symplectic manifold and let pr : RM7>M be 
the projection. For a given regular function H : Resse derine: Ene 


2-form 
re x 
(1) w= pr wt dH A dt 


Then (RxM, w) is a contact manifold ([1]). A motion of the system is 
a nonparametrized integral curve of a vector field satisfying X44 i S Oe 


If it happens that w = d6 then 


~ ~ 


w = dé 
(2) 


D 
I 


* 
phe lel Gan lah folje 


(d) Evolution space; space of motions. Denoted by P and M 
respectively, these were introduced by Souriau (see for example, [48]). 
Evolution space P isa generalization of the contact manifolds of (b); 
for a Hamiltonian expressed in homogeneous coordinates (time is included 


* x ; 
as one of the coordinates), P = T (RXQ) . P= R*XTQ is the evolution 


te . U J i 
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space of a (Hamiltonian) contact manifold. “ Defined in such an evolution 
space is a closed 2-form W of constant rank < dim P , that is, a pre- 
symplectic form. W defines a characteristic distribution p > ker Wp) < 
Mats which is integrable since dw = 0 , and thus p> ker w(p) isa 
foliation. Associated to this characteristic foliation is the quotient 
set M= P/ker w of leaves (that is, maximal integral submanifolds of 

p > ker w(p) ). Denote the projection P *M by MT - to each point 

p €P,M7 associates the (unique) leaf passing through p . By requir- 
ing that tm be both open and continuous, a unique topology on M cs 
defined. On the other hand, the existence of a nee ever ico le structure 
on M such that M becomes a quotient manifold and 7: P7*M a 
submersion, is not as easy. Call a submanifold U of P containing a 
point p€P transversal to the foliation ker w if at each point 

ah ewe Het is supplementary to ker (Dp) @sathat 1s, tee = ee ® 

ker W(p) . U is a section of ker w iff U meets each leaf of the 


foliation in at most one point (or not at all). 
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Then, the foliation ker w is said to be sectile (or regular) iff for 
each peéP there is a submanifold UcP such that peU, U isa 


section of ker w , and U is transversal to ker w. 


Theorem: With the above definitions, M = P/ker W is a quotient mani- 


fold iff ker Ww is regular. (For a proof, see Palais [37], page 19.) 


It the theorem holds, there exists a unique symplectic form w 
on M such that W = TO  tigteis ewe sonpere is the pull-back, of Ww 
on M. M is the space of motions; each point of M is a leaf of 
ker W , that is, a maximal connected integral manifold of the distribu- 
tion p> ker w(p) . In effect (P,w) is analogous to the example 
described in (b): WwW describes the dynamics on P 3; one solves for 
(local) vector fields Z on P satisfying Ziw=0 (iff Ze ker W), 


and the solution manifolds are the leaves of ker W . 


(M,w) then contains all the canonical structure of the system, 
for example, the Poisson bracket structure on the space of all observa— 
bles of M. An observable f on M is defined to be a real-valued 
(smooth) function on M , so we show how to define a Lie-algebra operation 
on c (M) . Suppose f,g € c (M) > recall that W defines an isomorphism 
oF : *(M) > ¥ (M) of vector fields to 1-forms: 62) Sn ll) x (™) 


for Ze #(M). Then 


-1 
(1) {f£,gh = (oe (df)) A deg 
is the Poisson bracket of f and g . Local expressions and further 
properties of Poisson brackets are found, for example, in [1]. Canonical 


structures of classical systems is of contemporary interest because of its 
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role in canonical quantization ([21], [10]). Although the theory of 
finite dimensional systems is very well developed as well as the main 
interest being field-theoretic, a thorough analysis of finite-dimensional 
systems is often indispensible in the study of infinite-dimensional sit— 
uations. The latter are symplectic manifolds modelled on (i.e., coordi- 
natized by) Banach spaces. A lot of the finite dimensional theory 
carries over unchanged, but care must be taken, for example, with the 
different topologies involved. For a recent survey, see Chernoff and 


Marsden [7]. 


84.5 Space-Time Model of a Free Rigid Body. 


(a) We now take the principal bundle of G - frames, P , over 
space-time V as evolution space for the rigid body (motation - 84.3). 
The canonical system (P,w) is defined in this section, and the calcula- 
tions of the canonical distribution ker w (84.4c) and the associated 
equations of motion in space-time are carried out in the following sec- 
tion (§4.6). The treatment of the relativistic system given here is 
essentially that of [28] while the Galilean, in this formalism, appears 


to be new. 


(b) Space-time. Let V denote a four-dimensional differentiable 
[ee] 
(C -) manifold. The following properties of a space-time are not unreas~ 


onable: 
(i) Topological: paracompact, Hausdorff, connected; 


(ii) Orientability: V is oriented (i.e., there is a globally 


defined volume element oO on V ); 
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(iii) (Time Orientability); there exists a tangent line bundle 


om Wo 


The assumptions (i) are not special; they appear to be standard 
in the literature. It may be pointed out that paracompactness implies 
the existence of, first, a partition of unity on V , and second, a 


(global) Riemannian metric on V_ ([59]). 


Property (iii) is equivalent to the global existence of both 
a Galilei and a Lorentz structure on V (cf. Steenrod [50], page ZOII% 
Even if (iii) is not assumed, it will hold for any compact manifold 
whose Euler characteristic is zero ([50], page 203) sor tor any non= 


compact manifold ([61], Theorem 4.8). For more details, see [60]. 


(ii) and (iii) together in the relativistic case imply clatcue S/ 
is space and time orientable. Physically, it appears that time-orienta- 
bility implies space-orientability. This is argued on the basis of 
elementary particle reaction experiments and also the CPT-theorem (see 


Hawking and Ellis [19] for a discussion and references). 


It may be noted that (iii) additionally is equivalent to the 


existence of a globally nonvanishing vector field on V 


In what follows V will denote space-time with additional 


structure; 


(iv) a first-order-flat Galilei structure with a Newtonian connec- 


eaayam le ey aye 


(iv)' a Lorentz structure with the induced Levi-Civita connection. 
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A free motion of a test particle will be taken as geodesic 
motion in both cases. Here comes an essential difference between Newton- 
ian and Lorentzian physics: The connection in general relativity describes 
only the gravitational force; in Newtonian space-time, the connection may 
include other force laws in addition to that of gravitation. Referring to 
4.3 (10), we may have a conservative force (gravitational; electrostatic) 

A AT 


It? palbe gam ey dno for some potential $ , or an electromagnetic-type 


. A = " * iA] 
force: Ve 2y raat] (for vadivector potential, Ar: 


(c) Space-time configuration of a rigid body. The space-time posi- 
tion of the centre of mass of a rigid body is a point xeV. As in 
§3.2, the orientation of the body is described by a triad of space-like 
orthonormal vectors (e,) . Orientability of V allows us to choose a 
unit time-like vector e_ at ae) In the Lorentz case, this is unique; 
however, the condition in the Galilei case for ey to be time-like unit, 
ae wy = 1 =, has a three-dimensional (affine) subspace of a at each 
x as solution. Uniqueness may be regained by requiring 2. to be the 
unit tangent vector to the worldline of the centre of mass, however this 
is somewhat unsatisfactory. Alternatively, a unit time-like vector field 
u on V_ may be assumed to exist and required to have the coordinates 
6 in an adapted coordinate system (4.3(9)); this assumption is actually 
built into the formulae 4.3(10) ([28]). Im any case, we assume a config- 
uration of the rigid body is represented by a point p = (x,e,) in the 


bundle of G - frames (84.3). 


It is clear that the choice of G(V) as configuration space 


for the rigid body is appropriate. Not only do the physical quantities 
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have natural interpretations in terms of the geometrical structure 
(explained more fully below), but it is easy to demand some degree of 


covariance of the theory. 


(d) TG(V) as evolution space. Recall (Appendix) that TG(V) is 
parallelizable;{e ,E,,F,} forms a global basis of xX(G(V)) while 
{6° 2 Fo} is the dual basis of x" (c(V)) . It was also pointed out 
that Ey generates so(3) at each point in V , so wn, is again 
angular velocity at xe V in the body frame. Denote coordinates on 


TG(V) by Gor ceies ay lg) 


~ 


To see what to do next, examine the canonical one-form 60 , 
discussed with respect to time-dependent canonical dynamics in 84.4b, 


for the particular case of a free rigid body (no external forces are 


present): 
tea A A.B A mb 
(1) § = 7 AMV yV + Tape ade my ,0 1,30 E 


(see 3.3(8) and 4.4(2)). We follow Kunzle [29] by taking as Cartan form 


Oo AGB 
(2) 6 = MO - Tape E 
where 
1 B 
3) M=m+5 1,00 ot Nasik ater 


in the relativistic case only (cf. 84.1). Here, m denotes the mass of 
the body, measured in a frame in which the body is at rest, both trans- 


lationally and rotationally. 
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As a motiyation for (2), we first remark that, by assumption, 
the model of rotational motion developed in 82.4 is valid at each space- 
time event; whence, we retain the last term of (1). To justify the 
translational part of § , we recall that if the two-form w is exact 
(as it is here), w=d0 , then associated to a dynamical symmetry of 
the system, generated by a vector field € (so that L,w=0), is a 


E 


constant of the motion, fe » given by fy = -§-4 6. Now, it is 
widely accepted that p = Me, is the four-momentum of an extended 
rigid body ([12], [13], [29]). p can be considered as a generator of 
space-time translations in the p - direction; under such a symmetry, 

2 


the corresponding constant of the motion is ee =ptoe “Me + g=—M , 


as would be expected. 


To complete the definition of the canonical system (T Lor(V),w) 


in the relativistic case, we define 
(4) w = dé 


where 6 is given in (3). The calculation is carried out using equations 


4.3 (18) (with v= 1 ) and results in 


ee) a,b Oa 
(5) we ZS Ry, OA +S, ast an ui A0, 
il Rae roan: ah gs 
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(6) é 


In this case, Sy af is the spin angular momentum in the body frame 


(eompares2a4 (24)) 4 


We return to the Newtonian case and attempt to mimic the above 
procedure. If we try to use (2) to define a Newtonian Cartan form 0 , 
we immediately run into insurmountable difficulties. For, not only does 
Meeceducastogm (4 constant), ¢so sdMs= 0.; but also, for the Galilei 
structure, d0° = 0. On passing to w=d06_, we lose all translational 


information, and the procedure therefore breaks down. 


What we will do then, is to postulate a presymplectic form 
in the Galilei case. To give a reasonable justification for our choice, 
we find it convenient to consider simultaneously a contact structure on 

3) : , 5; Lee pale 
E = RxT(R XS0(3)) which we coordinatize (globally) by (t;x rE ,5V AVEO 


and a flat Galilei structure on Ro which is defined by 


0 
Va oF 
C7) j 
' g@le O ABB 
Yo age O49 6, 
The Cartan form on E is 
il lf B iP oA 
(8) 9= 5 (mv, + eo )dt -mvpdx -S,E , 
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(9) W = my .dv' Adt + Si der ade - mdv,Adx” 
~B all ~A ~B 
~ Tin AQ AB sree Sy yee 3S 


where we have used 2.4 (18). Now, we are used to thinking of pt as 


ie 
the velocity of a free motion, that is, = a ; otherwise written, 


(10) OE Glee an 


The exterior derivative of (10) yields dy Ade =O, 850s thal thesrirst 


term of (9) disappears and we are left with 


(11) W = -m dv ,dxe FIBA at nat 
A xB il Ned 
“Tap UGE NO Ge > ap E AE 
Now consider Ro with the Galilei structure defined by (7) 
and flat connection, [ = 0. From 4.4 (18) we calculate that 
iN [ =A nee . : ; 
de. i I i I Nelabel ahghen cone eo as a transformation to the velocity 


along the motion, we make the following correspondence between E and 


Gal(R‘) : 
dx <r e) ee 
(12) ag ae e, a 
dt <> @° 


Using (12) in (11), obtain 
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oO ~A Oo 
(13) W = m 0 AF + S 4 dy Ad 


~A ib ah @B 
sts TAR E adsl + > Sap E AE 


The choice for wW in the general Galilean situation will be the simplest 
one reducing to w in the case of force-free motion and flat Galilei 


structure, as well as being consistent with the relativistic form (5). 


Hence, the model for a test rigid body in space-time is the 
canonical system (G(V),w) where w in both the Lorentz cnd Galilei cases 


being written as 


2 oy aap aN iN Ae) 
el S Runab AG) MO AF oe SA df} Ag 
SA Biase) ~ho~B CIN Geel 
+ IB E Ad{? + > Sap (E AE -VF AF ) 
(14) O , Galilean case 
Vv = 
ie eLOSeciCZmcace 
aL B 
Memtvs I, 1 


§4.6 Canonical Distribution on G(V) ; Motion in Space-Time. 


(a) We defined in the last section the canonical system 
(TG(V),w) and now proceed as outlined in 83.3 to calculate the distri- 


bution p > ker w(p) 


Let 226 Deu be any vector; say, 


,4 







aS a Pa < io al tle 


ik) Aes ad Elite anh rivet [uo ie oily fy . hod | oh 43, 


ae 
ind Jaks bes 9/0 oo agi Bq! ies oir eae aaa aa 


nwo?’ \£9é¢0vtinias ec7 safe gaeeeheaes saled oa Lis aa 


ak 
= 7 : 
amii-souge at ¢Ybow inl d Feng £308 fate eae aor 


iG a ‘ rt 
Bad $45 Ps * A 


70. 


(1) Z= vie + PE oe et a yA 


Oe 
. ov" au. 


where the coordinates on G(V) are (x7 2) and those on TG(V) = 
a 


G(V) x Ro xg are Gere. ea: ra) . By definition, 
ker w(p) = {Ze T 16 (V) Zs} O(p)e= Ol . Now, in terms of (1), 
7 eh gn 
(2) Ze gta-—CS Piers Stak nae 
1 N TA 
eh ae Runar ¥ 
B BL +A 
-(1, 5° + S ape Re + Qt, + Vv aN 22 )F 
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O A 
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where we have made use of the fact that fe ,E,, an ; ith t a 
a Ey Ea a of 

Ov re) 
and (92, EO Fos dv? ag,au} are dual orthonormal bases of BNI and 


k 
Ea NAC, respectively, at each p € TG(V) 


Therefore, Ze ker Wp) iff Z2w(p) = 0 , iff 


(3) pA = y° of 
(4) My = -V aa 
A o ~AB C 
0) View Sac i?) 
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(7) 0 = See 
where (5), implying 

A 
(8) | oy One=s Ole, 


is used to calculate (7). Note also that (4) implies 


(9) S.v = 0 


Substitute (6) into (4) to obtain 


A 1 MN reba 
(10) Mv - v 7 8 Rune BY Saaee= Oa). 
an equation for oe in terms of the coordinates and v° OnLy wil 
the Lorentz case, the solution of (10) is facilitated by first transfor- 
AB AB i 2 3 


ming to a frame such that S” =s 653 ESOL vAnges Olea eee Usene t\) oo 


and applying the inverse transformation. Equation (9) implies that vA 


has no component along Ss . We find then, that 
Ane ele ee? 71 Ol MN BA 
Gia) Vie Si eee U eS eR ee 
with 
a 129 MN AB 
(12) AN SESE as Z MM S RNABS 


The assumption that A not being zero places a restriction (in the 
Lorentz case) on the magnitude of the spin ([27]). Physically, the test 


particle concept breaks down if A= 0 (for fixed non-zero curvature) , 
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so that the condition A #0 is not significant. We will later consider 


(12) as a guide to approximation in Schwarzschild geometry. 


A wes aioe Lee LENG » whose components (1) are given by (3), 
(5), (10) and (6) with Ue : ba : yA remaining arbitrary, lies in 
ker w(p) ; whence dim(ker w(p)) = 8 and SN CL NUON Sees w(p)) = 12 
This holds at each p = (x ,eo5Vv%, U)e€ TG(V) (subject to (12); we 
interpret the evolution space now to be the open submanifold of TG(V) 
defined by A+0) , so that ker w is a foliation and the space of motions 
M = TG(V)/ker w is (assumed to be) a 12-dimensional symplectic manifold 


corresponding to the traditional 12-dimensional phase space of the rigid 


body. 


(b) Motion in space-time. The motions of the system are the maxi- 
mal integral submanifolds of the foliation p > ker w(p) . The equations 
describing the translational motion of the centre of mass of the body 
and the rotational motion are the equations of the projections of these 
integral submanifolds into space-time V , and G(V) . We begin with 


the following definitions: Let 
(13) A POA) = GA), e2A);v7), MO) UA) © TE) 


be an integral curve of a vector field Z in ker w. By definition, 


A + p(r) satisfies 


(14) ar pA) = ZA dp (A) 


We use the notations 
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The following space-time components of the corresponding body-frame 


quantities will be needed: 


(16) oP = ey e 
(17) ga ae 
(18) gl? en ef oo ee = cas e! 5° 
(19) 21 7 94 8, es = €vey5 ev a 


We are interested only in the worldline of the centre of 
mass, > x (A) , the corresponding development of the body frame, 
. QO 
Na en () , the angular velocity, A> QA) weandmthe spin, A->7S (A) 


Explicitly, then (13) is: 
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Buty trom: (20) ,5 50. s= OF v3; using this in (6) we get 
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(Note that Sees or = Seah ee , in view of (7).) Similarly, 
de® e 
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Note that (in the Lorentz case) ae e = oe se = Ons etaus 
A Go APa8 O d : 
a e = Sova and (22) implies ee = Q Substituting this 
A Qo A do A ao Aaa 


value for Qa into (25), and making use of (3) and (19), we find 
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However, noting that wh a. =] and y or = by evel? = 0 eeine the 


Galilei case, (26) can be equally well written as 
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is the projection operator onto the instantaneous rest space orthogonal 
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Substituting (24) into this equation and using 
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we find that 
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whence, recalling (22) and (28), 
=O 36 © 0s) Y 
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(30) 8 View a , 
where os = eves Tee (30) is the space-time form of Euler's equations. 


The last quantity we wish to consider is the spin ae 
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This essentially completes the description of the motion of a 


rigid test body in space-time. 


We close with some remarks and an example. 


(c) Remarks. 
(i) Due to the arbitrariness of the parameterization i + p(A) 
G13) v is arbitrary; if . is chosen to be proper time in the 


Lorentz case, or 


is determined by 


(ii) in the 


is automatically 


fact that (en) 


Oo 4 F 2 5 
A=t=x (absolute time) in the Galilei case, ve 


Ca 


a = 0 


Lorentz case, 
Oo 


the condition Sy or equivalently, 


one of the so-called spin supplementary conditions, and 


satisfied here because of the definition (17) and the 


B 


hye 0% 


: O 
is a Lorentz transformation (so that = sei 


From the beginning, these supplementary conditions were recog- 


nised as being necessary to give a complete set of equations of motion, 


along with those 


ique, see Dixon [12]). 


for linear velocity and spin. (For a historical crit- 


Specifying such a condition on the spin has the 


effect of determining the centre of mass of the body, which, in relativ- 


ity, isa rather 


tenuous concept. A comparison of two conditions, 


s?° = 0 and gob Vas O , made in a practical setting in [4], shows that 
there may be a significant difference in the resulting equations of 
motion. 
AgLie 
(iad) the quantities Mi= mit Vv 58, ci and s = (SS ) / = 
Gina's are conserved along the solution curve. This follows by covar- 
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iantly differentiating M and S , and using the equations of motion, 


(30) sand (31). 


(d) <A simple example in the Galilei case can be obtained by consid- 


ering a flat Galilei structure with the induced metric oo, = On, . The 


condition that a connection [ on Gal(V) be Newtonian, dk = 0 


(cf. 4.3 (13), is equivalent to the local existence of a 1-form 


LL 


= oO Soc tes = i 
A Ay dx such that kK 5 dam ob eae Ja Define 
HP = ee 0 AQ and gt = 0 A, .| Reterring to 4.2)(10), we see that — 
the nonvanishing components of the connection are coe a -E) and 
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(e) Fermi transport of spin. Now that the equations for the 
centre of mass worldline and spin are known it is useful to see how (31), 
the spin-precession equation, can be interpreted. First, we recall the 
definitions of Fermi-Walker derivative and transport. Let u=u oy be 


: : a, : 
the unit tangent vector to a timelike curve A> x (A) in V and 


let €= ae be any vector field. Then the Fermi-Walker derivative 





Of = 1's 
ee Wee We Gal 
(35) CHa ‘i. 2 
: ae + us" a oe oe Lor 


—E is said to be Fermi-Walker transported along the curve iff 64e = 0 
on the curve. Fermi-Walker transport preserves scalar products, so that 
an orthonormal frame at one point remains orthonormal when Fermi-Walker 
transported along a curve - such a frame-field is called a Fermi frame. 
A Fermi frame is the natural space-time analogue of a non-rotating 
frame. If (e,) isva Fermi trame with ec. = Us | then (e,)) is a non— 


0 A 


rotating triad in the instantaneous rest space orthogonal to u. 


Applying (35) to the spin vector Ss” we get 


f 
| 0 » Gal 
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(36) Cops) ie ie . ; 
i. ge “2 usu" = ue a Lor 


where (31) was used for the Galilei case. We obtain immediately the 
result that spin is Fermi-Walker transported in the Galilei case. To 

f a, 
compute (36) for the Lorentz case it is necessary to first obtain eo 


in terms of a in order to make use of (31). This is done by inverting 


equation (21). 
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All considerations now refer to the Lorentz situation only. 


Define 

(37) u = ar , and 
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(38) oe Smt 


Thessolution of (21))) for ae in terms of ue 
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decomposing €y as 


Q SO a S0u 


(39) e. =u’ +e ,; gu = Oe 


Substitution of (39) into (21) yields 
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from which it can be derived that 


(41) a= s0no a= 0. 


can be obtained by first 
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Thus, in frame components (40) becomes a 3-dimensional problem: 
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which can be solved using the trick employed to derive (21): 
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Returning to (39), we see that 

ee es Re ee) m9) 
(44) . en =u 5 M S Beat 


Now, define 


epelnt oY DEO 4.08 
(45) A" = a M 


| 
ae 






5 hee | 
ay bait aw “s 


: | 7 - 
| | at) 
tmis biviwel ed dee ss loki on : ; 


Sm 2 “penogmn. sma’ ak 7 
} & | eae) See 


ar 
; ¥ 
- gt te a 2 ae aie ¥ 4, 
ee in 





80. 


Then using the definition (37) of u together with the result (44), 
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(46) WI 1-9’, 
VV 
a 
from which it follows 
22200 a -2 UV 0,00, 
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(47) (Gi ) @ey =u *M “S Root S 


(We shall see that in reasonable situations A will be small compared 


to c= speed of light = 1.) 


Returning to the case in hand, note that (21) and Se 5 = 0 
imply 
a. 
(48) S ee 0, 
so that 
(49) See eh 
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Equations (36), (48) and (49) show that 

(50) (88)° = (69 + ud, )e 

For brevity, introduce 

esis) Ree ee ouk (ees 

so that (31) becomes, using (47) and (50), 
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Thus, in the Lorentz case, spin is not Fermi-transported along the centre 
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of mass worldline. In 85.9 we will see in a typical situation how 
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CHAPTER V 


Constrained Motion of a Rigid Body 


G5. sintroduction. 


(a) Before proceeding with the main subject matter of this 
chapter, we make a few comments of a general nature. The usual discus- 
sions of constraints in mechanics often don't make quite clear what 
issues are involved and how. We believe that a discussion based on a 
geometrical formulation of Lagrangian dynamics can provide both arnacur— 
al conceptual framework as well as the necessary apparatus for dealing 


with concrete problems. 


Standard usage of the verb "constain'" implies that something 
is being constrained to something else; otherwise said VConstrain 15 
an operator. In our context, to a given dynamical system, this opera- 
tor assigns a new dynamical system in a manner which depends on the 


particular description of such systems. 


For example, the evolution space of the traditional Lagrangian 
mechanics is RX TQ, Q being the n - dimensional configuration 
space of some dynamical system. An associated constrained dynamical 
system is then understood to be obtained from a sub-fibre Bundle WN _ over 
R of the fibre bundle Rx TQO*>R: For each t € R we assume that 
there is given an (m< 2n) - dimensional submanifold N, Offs sO so 


that N is obtained as N=u Ny . Reducing the evolution space from 
IE 


Rexe iQ, tow Nowwillshavegan effect on the dynamics, in that constraint 


forces, a priori unknown in general, will now appear in the equations of 
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motion. These constraint forces are often thought of as being the 
agency which restricts the motion to N. If N is a product, 

N = RXM where M is a sub-bundle of TQ, the constraint is said to be 
scleronomous (time-independent); otherwise, it is rheonomous (time- 
dependent). Furthermore, if each N. ¢ TQ is derivable from constraints 
in Q, N, = TC, , C, © Q being a submanifold of Q, then the con- 
straint N is called holonomic; otherwise, it is called non-holonomic. 
The generalization of this example to a homogeneous description of 
mechanics, the evolution space being T(R*xQ), is straightforward, and 
we consider it to be done with the understanding that Q now contains 
the time coordinate. The constraint manifold N is then a submanifold 


ene AMO), 


The case that dim N < dim TQ (holonomic or velocity-dependent 
. constraints) is often described in terms of a smooth map 9 : TQ > a 


where m <n = dimQ; that is, N= enero). LP oO, 


ee 
‘Deeps 
area ie is surjective for all peN, then N is an embedded submani- 


fold of TQ and has dimension 2n-m ([59]). 


When N is holomonic so that N = TM for some submanifold 
McQ, describe this M (since all questions studied here are local 
we will denote again by M a certain open submanifold UcM) bya 
submersion ¢ : Q> R" (an), M = Soe, say. A vector v é me 
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Leeisectearathat the map 0 = (04,6 )° =: TO°> R" is a submersion 


whenever is (compute the Jacobian). 


(b) Historical development of constraints. An important category 
of constraints are those arising from systems described by degenerate 
Lagrangians. We now try to indicate the difficulties involved and what 


they lead to. 


The Hamiltonian formalism of mechanics appears to be the most 
convenient as a point of departure in the quantization of classical sys- 
tems. As well as being the traditional way of proceeding from classical 
to quantum mechanics, the Honiton tan formalism gives us a canonically 
defined Lie algebra of observables (cf. 84.4d) which is the basis of 
some modern prescriptions for quantization (for example, Hermann [21]). 
On the other hand, while the requirement that a theory be relativistic 
is easy enough to formulate in terms of a Lagrangian and its associated 
action principle, this is usually not the case in a Hamiltonian formula- 
tion (this becomes especially apparent in relativistic field theory where 
the Hamiltonian formalism treats the time coordinate and space coordi- 


nates very differently ([11]; [42] page 287)). 


The Mroceatre outlined in 82.5 which constructs a canonical 
Mioomeian structure on TQ depended on the Lagrangian L being regu- 
tape (Zoo Clo))) to pull "back the canonical ‘syiplectic orm Ww, on T 0 
Come Oi [Lele Leeshowimuiac di “the tipperdcetivativew fin asa 


diffeomorphism, then the Hamiltonian and Lagrangian formulations of 
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mechanics are equivalent - the projections of the solution curves, from 

* ° e 
IQ and TQ respectively, into Q coincide. In this case, we can 
readily combine the advantages of the possiblities of relativistic covar- 
dance and canonical quantization by taking the Lagrangian formulation as 


basic. 


Unfortunately, there are very few physical systems of interest 
10) 
for which L is regular. This, together with the need to quantization 
led a number of people into considering methods for dealing with such 


degenerate Lagrangian systems ([10], [11], [26], [31], [45], [46]). 


(c) The purpose of this chapter is to investigate the constrained 
motion of a rigid test body in curved (Newtonian or Lorentzian) space- 
time. The next section shows how to geometrically view the traditional 
formulation of Lagrangian constraints and leads to an inhomogeneous for- 
mulation of symplectic constraints. §5.3 outlines a homogeneous formula- 
tion of symplectic constraints while 85.4 applies this to the rigid body 
problem; unfortunately, the relativistic equations of motion are not 


completely determined in this case. 


The reformulation of the problem as a second order equation in 
§5.5 in the spirit of [29] (Appendix) and [1] not only provides us with a 
deterministic set of equations of motion, but also a plausible reason 
for the failure in §5.4. §5.7 examines some simple examples (Galilei 
space-time, Minkowski space and a DeSitter universe) of motion constrained 
to a worldline. Finally, after showing how to obtain the general solu- 
tion, we give a "practical" example: motion constrained to circular 


orbits in Schwarzschild geometry. 
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85.2 Geometrization of the Lagrange Multiplier Formalism. 


(a) In this section we assume all constraints to be holonomic. 
Although holonomic constraints can be eliminated right away by making 
coordinate transformations in the state space, it is instructive to use 
the geometrical interpretation of the traditional Lagrange procedure to 


do thas. 


(b) Let Q denote the configuration space of a dynamical 
system with dim Q=n. Following the notation of 82.5, TQ is 
the state space and L: TQ >R_ the Lagrangian (assumed to be 
regular). The motion of the system is assumed to be subject to the 


constraint defined by f : Q> ne (f is assumed to be a submersion): 


(1) ey oe yA = eee, =< 


A motion of the system is represented by a curve R>QO _, hex (A) 
0, 
Sucimcna ted com) ice tOm alo est A (qr? (A) : “4-()) , satisfies the Euler- 


Lagrange equations: 
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(3) vi(ry = A 
where the Ay are the Lagrange multipliers, functions on TQ . Equations 


(2) and (3) are subject to the constraints (1), as well as the compatibil- 
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These are essentially equivalent, in the holonomic case, to (1), but 


we find it convenient to state them both here. 


Now, being a lift, A> (a°(),v Q)) is an integral curve 


Omeay vector f2eld 9 Z° on elO0) asuchi that 
C5) z dq” = vy" 


fice TO. eR ei Ghanvecunct1 On voles (d.v7))) then 


dF 


(6) Z1d¥ = Z(F) = = 


by definition of ) > (ants) being an integral curve of Z. Using 


(6), (2) becomes 


A 
(7) ze) = 4 SE, 
OV oq 0q 
so that 
(8) z(Sbydg’ = ag? +d, af, 
yt 9 oy A 
q 
* 
since ae = oF dq” Define the 1l-form 6 € ¥ (TQ) by 
oq 
(9) ee oe. dae 
OV 
(Ch; 52.5) and also note that 
OL 
(10) L, (0) = 2(=)dq™ + L, (da) 5 
hv OV 


where Ly denotes the Lie derivative. Use the definition (9) and the 


formula (10) in (8) to get: 
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We recall a standard formula: if a is a k-form, 

Gi?) L, (a) = Z-4 da + d(Z4a) 

Applied to (5), 

(13) L, (dq) = Z 4 d(dq*) + d(Z4dq") = dv 


The 2-form w is defined by 


(14) w = -dé : 

hence, 

(15) L,(8) = -24 w+ dv a 
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is the energy (cf. 82.5). 
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If we denote by M the constraint manifold in Q , 
A 
(18) Meese (qs £(0)F="07 ; 
the "lifted" constraint manifold N in TQ is just 
A =A 
(19) Nese Meee (Game Osc at (dq) o— Omega (as) = 0}; 
where we have defined the map f: TQ > Ra by 


 ( = ip) aes nuh tebe 


(20) 


The condition that the motion stays in N is that Z be tangent to 


and is equivalent to 


N 
be 
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Fh 
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Let 1: N+ TQ be the embedding of N as a submanifold. 
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. * * 
(22) Lido) = di a) 


We can State our results on N intrinsically as follows: From (16) 
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where W = 1 W (ei eZ eo C19) 
(c) Now, let us suppose that we are given the system (23) on N 


defined by (19). 
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(d) We can summarize the above discussion by the 


Theorem: Let Q bea configuration space dimQ=n ,TQ state 


space and L: TQ>R a regular Lagrangian describing a physical system 


subject to the constraints e Reg Av—sls wasn pwaeres, £ < 0 > R™ 
is assumed to be nondegenerate on M = {qe QQ: f(q) = O} . Let 
N= T™ . The following two descriptions of constrained motion of the 


system are equivalent: 


(Pee the motion 1sea curve: 9). > (q7 (A) vo (A)) in TQ satisfying 


the Euler-Lagrange equations 


ey Soe 
dir Jy sae A oe 
ay = $© @) 
dA 
subject to ea) = 0 


where the rn are functions on TQ to be determined. 


The dy are called the generalized constraint forces. 
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(i41)") the motion is’ an integral curve of a vector field Z on N 


satisfying 
ZAVD —“dke= 0 (on N) 


ne Aa * 
Where si anda teware associated tor tl “as in §2.5. 


§5.3 Geometrical Theory of Constraints. 


From now on, the mechanical systems we will be dealing with 
will be described more generally. By this we mean that no longer 
will a Lagrangian describe the dynamics; but, we will assume that a 1l- 
form 9 such that w= d0 is a (pre-)symplectic form, or even a given 


presymplectic form w , is the basis of dynamics. (Cf. 84.5.) 


Constrained dynamical systems will first be described abstract- 
ly; then we specialize to the case that constraints are holonomic. It 
may have been noticed that in the last section the geometrical formalism 
induced by the Lagrangian could not describe nonholonomic constraints. 

A geometrical theory of mechanics which appears to be especially suited 
to handling nonholonomic constraints has been developed by Vershik and 
Faddeev [58]. Essentially they use a geometrical structure on TQ it- 
self (where Q is a configuration space) consisting of a mixed-rank 
second-order tensor t on TQ , a 2-form 0 associated to a Lagrangian 
PeeOnweLOmavia At . andr aavectOneiaeld x on TQ defined by means of 
Q and L , whose integral curves represents the motion. (Klein [23] 


also considers this particular geometrization of mechanics in a more 


general setting - he recognizes the geometry as being an almost-tangent 
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structure). Constraints are then introduced as a submanifold M ¢ TQ 


and related to the mechanics by T and the sub-bundle TMc TTQ . 


Following 84.4d, we define a dynamical eyeren es be a pair 
(P,w) where P is an n - dimensional manifold, the evolution space of 
the system, and wW is a presymplectic form defined on P. A pair 
((P,w) , 1: P' > P) consisting of a dynamical system (P,w) together 
with an embedded submanifold P' , 1: P'’ + P being the embedding, is 
called a constrained dynamical system. Unless explicitly stated other- 


wise, P' is assumed to be contained in P as a subset. 


Given. a constrained dynamical system ((P,w) , 1: P' > P) 


we can define the following four vector spaces at each p € P 


Ca) we ‘}the tangent space to P at p 
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Using the dimensions of these vector spaces (assuming dim P =n _ fixed) 


we can define two integer-valued functions on P 
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p > dim N -— dim K 
Pp Pp 
We assume that these functions are constant on P (an indication of 


what happens otherwise is found in [45] and [46]). We will denote by 


K both the sub-bundle U Ke of TP' and also the distribution 
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x 
p! > K ,.. Since “dim’K jis constant on Bs @keeckern (1.9) esisea 


completely integrable distribution and hence defines a 


P' (see 84.4d). Points of the quotient set 


(or histories) of the system. 
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More concretely, let Q bea n - dimensional (extended) 
configuration space of some system and let P = TQ . Suppose that 
‘fe ORS R™ (mn) is regular on C = fei) Sh cece Oly tq a= 0} 
Denote the component functions of f by = Osos RAIS. Lt). mee And 
Vetern(g bee =)0,... ,n-l)ee be a-local coordinate system at some point in 
Q. Then P' = TC is described as the set of zeros of. the map 
fe eae Oa R xR” : £(q,v) = ap) eon =v- EDD OLE Ocal ys 


—A Q ye 
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it > e 
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oquee 


The conditions (4) can now be stated as 


a A —A 
=) = 
v~ w Ay df + Uy df 


<p eee Ss i) eee = 


(5) 


aN and Uy to be determined 


85.4 Application to Constrained Rigid Body Motion. 


(a) The general discussion of the last section and in particular 
the equations 5.3(5) will now be specialized to the evolution space 
P = G(V) and the presymplectic form discussed in 84.5 and 34.6. We are 
interested in seeing the effect of constraining a rigid test body at its 
centre of mass to a specified worldline and in particular a geodesic 


(with respect to the connection on G(V) - see 84.5). It is expected 
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that no effect is seen in the Galilei case, but that a precession of the 
spin axis will occur in some (general) relativistic situations. Unfor- 
tunately, we don't see this now but have to wait until 85.8 and §5.9 
since the present formalism leads to equations which do not completely 


determine the motion of the system. 


(b) The definitions and conventions used here are those found in 


84.5. A qersiiblar represented by )\ > wv) Poca ily a teeeVelecal 
be locally described by {xeV: meres = 0} where for example 
i Ik i 
( o (x) = x = W (g(x°)) 
(1) 


xo = P(A) «iff A = B(x’) 


4 ip 
We denote the lift ae of oe into G(V))*again by > (so, 


6’ (x,e,) = > (x) ). 


The equations 5.3(5) are now written as: 


(2) = 0 

(3) oe = yt Ve be = 0 
(4) - 74 OS ali = 0) 
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The vectorfield Z has the form §4.6(1): 


z -v'c t+ PE +e, +e + oh Sey S 
av" 9M aU 
Introduce the notation: 
r_oT ey aoe 
(EO F = 6 =e ap 
a la a 
Ox 
(8) e = al 
Oo 
he OP plea! 
(9) (Gn) = ©) 
(10) af -cr fa 


We have assumed that $ is regular, hence (F.) has maximal rank, 


Denote by M the submanifold of TG(V) defined by (2) and 
(3). Using the formulae 84.3(17), §4.3(18) we calculate the differentials 


of the constraints: 


(11) pues F BE 
mee ase S if =B 
(12) dd =v FY gy = Fy Vv, €& RB E 
Owns Ralls age 
+ (v F. + Vv £ Vp )F + F dv 


Equations (3) and (4) are computed to be F, v' =0 or 


yh aia, he , and F vy? =0 or vA = Be , respectively, and therefore 


imply 
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O 
(13) v4 = a 5 
O 
V 
while (5) works out to be 
Leper: it iy B i KL B 
(14) F Eo + (FR +) £ V,)Q - F ee P 
ig? ago a We oe 
ris 


Finally, Z-+4wW was computed in §4.6(2), so that the use of 


(11) and (12) enables us to express (6) as follows: 


(15) so" + 5 Spite ve = dp oS Up Flo v* 
(16) : IR y= MQ, = Ap Fy he is v" 
(17) 43? ah ha Ser Fe Gy, al 

(18) my + Vv so = up(vF, + \ fv) 
(19) Teen Een)” Sip 7 0 

(20) O= uF, 


In the Lorentz case (v=1) (13) to (20) do not provide enough 
conditions to fully determine the motion (the difficulty bcing that the 
A Acar: ’ Xe) 

component of Q along S is not determinable; the equation for ee 


depends on of » hence e, cannot be determined). In the Galilei case 


(v=0) we can obtain the equations of motion as follows: 
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Equation (20) implies that 
(21) Up =O; 


since (F.) has maximal rank. This implies that: 


(22) Me 
(23) f° BAO 
(24) ce ay ae ee 


Let A be a parameter of an integral curve of Z (cf. 84.6b; in parti- 


Clare 4,oCl3)eand 4.0014) ). 
(25) Vee ee dk Se ae 


determines en along the a priori given constraint worldline Ne OD 
; f A A Mon AL 

(Note that this determines Q :e = en QO 7 elving 0 in terms of 

(25).) The remainder of the equations of motion are obtained in exactly 

the same way as in §4.6b. In essence, constraining the centre of mass 

of a rigid body in Newtonian theory does not affect rotational motion of 


the body. 


The discussion in the next section will both give us a workable 
model in the Lorentz case and explain why this case doesn't work here, 


whereas the Galilei does. 
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§5.5 Reformulation of the Problem as a Second-Order System. 


(a) Many dynamical systems are describable as second-order equa- 
tions (SOE) on some manifold M ; indeed, the classical Lagrange equations 
were formulated as such in §2.5. The remainder of our discussion of the 
rigid body is based on the consequences of assuming such a second-order 
principle. ‘It turns out that the condition that the vector field Z 
found in 84.6a be a SOE imposes new Bone trainee in the Lorentz case, but 
results in no new conditions in the Galilei case. We begin with a dis- 


cussion of the SOE formulation outlined in the Appendix to [29]. 


(b) The setting now is the same as 84.5 and 84.6a. We assume that 
a motion of the system is determined by a second-order equation Z on 


To (Vee thatels a eZe= Seots the torm 
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(1) z=u(vee +0°E +UE,) +8 


where U , ED “ of ‘ ys are arbitrary (at this stage). We still require 


Chateecee eke We, 
(2) Ze w=) 0 : 


but substitution of (1) into (2) and using the expression 4.5(14) for w 
shows that (2) can be solved only on the submanifold Nc TG(V) defined 


as the set of zeros of the functions: 
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5) tite a ee 


where we have introduced the shorthand notation 


(6) R 


N]FR 


MN 
: Runab 


ab : 


Differentation of (3) to (5) shows that ca 5; r are independent on N. 


One additional condition arising from (2), 
(7) We ROH a0 | 5 


is a consequence of multiplying (5) by v, and noting that (4) implies 


A 


(8) as uS = 0 Son Ne 


Where TG(V) is a 20-dimensional manifold, the new evolution 
manifold N is 13-dimensional. The condition that Z be tangent to 
No, 2 = dy* = 0 and Zz 4 ar = 0 , allows the components ee ; ys to 


be determined up to the factor w and a Ts esti eciven in s4.0 (5). 


The pair (Reto) is seen to be a contact manifold (in the 
sense of §4.4c) since p+ Z(p) is the 1-dimensional characteristic 
distribution on N defined by vw . According to [26] the unparamet- 
erized integral curves of Z represented by > p(A) coincide with 
Piretenves OEmone LOLlataOnmepe 24D) Or rere ne on N. The curves 
AX > 1(p(A)) in TG(V) are contained in the foliation defined by 
ker w on TG(V) and the two quotient manifolds TG(V)/ker w and 
N/ker (1 w) are diffeomorphic (we are assuming in the first place that 


these quotient sets are manifolds). 
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However, our interest is in the submanifold C ¢ TG(V) which 
is the intersection of the constraint manifolds M (5.4(2) and 5.4(3)), 
and N ((3), (4), (5)). The next section carries out the constraint 


procedure outlined in 85.3 for C = MnAN ¢ TG(V) 


§5.6 Calculation of the Constraint Equations. 


The constraint formalism defined by 5.3(5) becomes, in the 


context of 85.5: 


(1) ro = 0 

(2) Ano 

(3) T=0 

(4) $7 = 0 

(5) § = 0 

(6) Z-2dr° = 0 
(7) Ze de eo 
(8) 7aaTs wee 
(9) Zed 70 
(10) Zo ab = 70 
and 


e. O A A r =I 
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Equations (1) to (5) we call the constraint equations, (6) tor (10) =the 
compatibility conditions, and (11) the foliation equations. The task of 
this section is to evaluate the constraint equations, to calculate (6) 
to (10) subject to the constraints and to write out (11). A few special 
solutions are worked out in the following section while the general 


solution is presented in §5.8. 


Explicitly the constraint equations are: 


(12) veal 

(13) Mv, + V Sy, tee 0 
(14) MU, - ve Ry = 0 
(15) o =0 

(16) ma - oi = 0 


Calculation of the differentials ar? : Ae ; at : a - 


aon is straightforward although tedious. The results are: 
Gey che coke 

(18) dE, = Mav, + V(v gS, + egg Ue 1p)dv? + 8, aU” 
(19) hy oS . SW oe Banh aN e 
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(21) ieee oh Fe je° ore et, Be 
s (F, a) i"v,)F + F dv* 
In equation (19) the shorthand notation: 
(22) ey (Geo taney) 
AB AB B oNrA oA B 
(23) aR = San 2 Rees eee oe UK = 
is used. 
Le 
Z = ve, + PnE, - OnE, a — +00 “s+ ys —. 
Ov ele oU 


denotes a generic vector field on TG(V) , then use of equations (17) 


to (21) gives for the compatibility equations (6) to (10): 


(24) Eo = 0 

(25) ME, + v(v,S, + Mula )o +v $5 ¥ = 0 
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(29) 


(30) 
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be Oe 


We use the expression 4.6(2) for Z-—-4+ w together with the 


differentials (17) through (21) to obtain from the foliation equations 
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if 
The notation oe Fy Op » Uy 


has been used in (32) through (39). 


ee eo 
= Fy Up : By = Fa Br and Die 


At this point, the necessary equations have their explicit 
form: we must solve (32) through (39) in conjunction with (24) to (28) 


ands (12)° to (16)% 


Before concluding this section we derive a number of relations 


for the Lorentz case (v=1). Equation (39) implies that 


g8 = -us? s*Pu, + s7(8,8")S" 


where 


(40) ea = § sh = : S cae 


(41) ch ee 6 


Multiplying S into (35) together with the use of the definition (22) 


A 
and equations (38) and (41) result in 


A — 
Therefore, 
(43) cae. “ie eee Ly Ge =a) 


It is a straightforward computation that 
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(45) a= 0 


§5.7 Special Cases. 


(a) V is a Newtonian spacetime; we take the equations derived in 


the last section and set w= 0 to obtain: 


Constraint equations 


(1) vo=l 
(2) vi = 0 
(3) MU Ree 20 
(4) g! = 0 
(5) fl = 0 


Compatibility equations 
(6) E> =F0 
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Foliation equations 


(9) Ba i Ba = On a 0 =O 
(10) Pie oo 
A _ _.o AB C 
Ps) 0) er al Sac" 
AgeecO ee A, ATE? 
C2) A =) (Ro at Gr oleae 
The main feature to notice here is that since a = 0 = Ba = 
Uy > there are no forces constraining the system to N. This essentially 


means that the condition that Z be a second-order equation is automati- 


cally satisfied in the Galilei case. 


The forces constraining the motion to M are given by (12) 


which is the same as 5.4(16) if one remembers the definition 5.5(6). 


NL ons a _ dx” a 
Since V =vv =0 , we again have v = as Z4dx” = 
Wee (5.4(25)9.9 That this is compatible with er = en of (obtained 
der a 
from St z- de. ) can be seen by explicitly calculating (8) in terms 


of the definition of at » 5.4(1), and then comparing with the deriva- 
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(b) V is Minkowski space; the conditions v=1 and R 


R 4 a. 
lsat it are used ‘in 85.6 to obtain: 


Constraint equations 
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(14) vin = 
(15) ut = 0 
cas et 
(17) fh = 0 


Compatibility equations 


(18) E72 = 0 
(19) U ek0 

A One Ae 
(20) Oa =F-\) Gr Palle 


Foliation equations 


(22) B, = “MY oy 

(23) Uy = aie Ty = eho See: F, : 
(24) Sale 

(25) ye Ce sere Shit 


o} was shown to be zero at the end of §5.6, while for we find 
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Unlike for the Galilei situation, the constraint forces 


A , : 
ae are non-zero in general in the Lorentz case. The formalism 
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of §5.4 does not take into account the conditions that the motions of 
the system satisfy a second-order principle; hence, the equations derived 


in that case did not fully determine the motion. 


The integral curves follow as in the Galilei situation; the 


equations for Q* and §° are the same as in §4.6 while those for o 


*Q we : 
and eo are given a priori by the constraints: 
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(c) V is DeSitter space time (or, a space of constant curvature). 


If kK denotes the constant curvature ([53]), then 
(27) RuBys = Kg 86 7, B68 ay? 
and K is related to the Ricci scalar by 
(28) Ke=— —— 8 


The components of the curvature tensor in a frame (e,) are 
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Recall (§5.5(6)) that 
(32) Rig = 5 ene ea" Ou NE 
Using this in 85.6(14) gives 
(33) U, = -M K 8,_v 
substitution of (33) into 5.6(13) yielde the condition 


(34) (m2 +K sv, ma) 
which we rewrite as 
(35) be: uw? s°K)v, =o 


At this time recall the discussion following 4.6(12); equation 4.6(12) 


evaluated in the current context is precisely 
(36) A=1+M “s’ K 
which is assumed to be non-zero, so that 

G37) vo = 0 


Although not evident at this point, we have come back to the 
same situation as in Minkowski space (example (b)). This is 
related to the fact that although V is curved, the curvature is 


completely isotropic and we expect no net gravitational force to act 
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on the body. It may also be noted that K = O is just Minkowski space, 


example (b). 


We take up the calculations again by examining the compatibil- 


ity equations §5.6(25, 26, 28). In this case; they are: 


(38) me* fn sty =O 
(39) KS Me 20 
(40) eg ti to 


Substitution of (39) and (40) into (38), and noting that (38) implies 


SA& = 0 , we get 


(m2 + PaO) = Or , 


hence 

(41) qY = =g* = 9 co ee 
Next, because v = 0 , we get: 

(42) a = 1), 

(43) ieee 

and it turns out that 
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(45) SAR 


With these results, it is easy to show that again 


(46) eek 
and 

A OC ~AB G 
(47) 6 =-y I Sac Q 


(d) Summary. We have examined in three examples of space-times 
the motion of a rigid body whose centre of mass is constrained to a 
given worldline. 
In all of these cases, constraining the motion of the centre 
of mass to the worldline has no effect on the rotational motion of 
the body. This is because equations (11), (25) and (47) which determine 
o” for each example are exactly the same as 4.6(5) and results in the 


unconstrained equations of motion 4.6(30) for Q and 4.6(31) for S: 


=a 38 Oe=a8 Y 
6 Qo =- 2 
B are Sey 
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Again, both M and S_ are conserved. (see Remark (iii) in 


§4.6c.) 


The fundamental difference between the Newtonian example and 
the relativistic examples is the appearance in the latter of constraint 


forces associated to the submanifold N (85.5b), but not in the former. 
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§5.8 General Solution of the Relativistic Constraint Problem. 


We return to the general setting of 85.6 and set v=1 to 
examine the general Lorentz case. A glance at the foliation equations 
and in particular 5.6(38) and 5.6(39) shows that most components of the 
vector field Z and the constraints a , By » etc. can be expressed in 
terms of Un s Accordingly, in this section we algebraically derive a 
system of equations which determine Uy in terms of the coordinates on 
T Lor(V) and the constraints al . The derivation is horrendously 
complicated in detail so we present only the highlights of the calcula- 


tions. 


In addition to equations 5.6(12) to 5.6(16) and 5.6(22) 


through to 5.6(45), we need the facts: 


(1) We 0) 


(2) S 


The latter equation shows that the problem is essentially two dimensional - 


+ BS ap . The follow- 


we can express HU, as, for example, Bye = ov, AB 


ing abbreviations are used throughout: 


(3) er She 
(4) y = SU" 
(5) n=vU, 9 
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The substitution of 5.6(43) into 5.6(38) gives 
(6) Bi Cy as 
A AB 


Equation 5.6(36) yields after some calculation 
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Substitution of (6), (7) and 5.6(23) into 5.6(34) eventually yields 
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If we now multiply 5.6(35) by Sap 
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and (6) being the major ones we get (eventually) 
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Now we focus our attention on the compatibility equations 
5.6(25) to 5.6(28). The idea is essentially the same as that used to 
derive 5.7(41): equations 5.6(28) and 5.6(26) are substituted into 


5.6(25) and from this we derive the desired relation for Uy, - 


First, though, we simplify 5.6(28) by making appropriate 


substitutions and reductions. It turns out that 
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At this point, all the intermediate quantities of interest, for example 


pa ; ae , have been expressed in terms of Un + After an extremely lab- 


orious computation, the equation for can be written in the form: 
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By contracting (20) with ve shoe Vp we obtain the system of 


equations: 
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Aa+BB=M 
(24) 
Cao+tDB=N 
where we have put: 
A = ee Tax al 
mes KU 
B=E=v Tax S Vy 
west K 
Cc =5 Ve Tax Vv 
(22) 
SAG KU 
D = Vr TAK S Vy 
M = san i 
— A i 
N=L SAT Vv 
Then: 
—1] 
(23) oa = A ~ (MD-NB) 
—1 
(24) B= A (AN=CM) 
with 
(25) A = AD - BC 


The expressions for the quantities (22) follow next. 
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Finally, note that 
(32) Sy L = 0 


It is not practical to write out (23), (24) and (25) in the 
general situations. What we have in mind is Cha terorsaspantacularepLoy— 
lem (22) is first calculated, then (25) and finaldya(23) and (24)3916) 


then gives and the other quantities such as ot (11) will folitow. 


Ma 
We therefore have an explicit algorithm for calculating the 

constrained motion of a.rigid test body in a general curved space-time. 

An example of this is given in the next section; however, we will only 


be calculating an approximate solution. 
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§5.9 Motion Along a Schwarzschild Circular Orbit. 


(a) The examples presented in 85.7 all lead to the result that 
constraining to a worldline (for those special cases) does not introduce 
any new features of the motion, namely, a precession of the spin 
axis. This was expected for a Newtonian or Minkowskian space-time and 
could be reasonably justified for a DeSitter universe. We turn now to 
consider an example in a Schwarzschild space-time with the hope that 
there is enough anisotropy present to produce some sort of effect 
(a Schwarzschild space-time admits a four-parameter group of isometries 


compared to a ten-dimensional group for DeSitter and Minkowski spaces). 


The example is that of a rigid (test) body constrained at its 


centre of mass to move along a circular orbit about the source producing 


the Schwarzschild field; one may think, for example, of a gyroscope 
orbiting the earth. We first give some details about such orbits then 
find an approximate solution of the equations of 85.8 and finally calcu- 


late the spin precession. 


(b) Schwarzschild orbits. We shall employ Schwarzschild coordi- 
nates ([35]; or curvature coordinates [53]) throughout. In this coordi- 


nate system, (t,r,@,6) , the Lorentz metric takes the form 


2r 2G 


(1) g = diag(-[1 - <2] , [1- 2)", x? , x sin’) 


where 2r 4 is the gravitational radius corresponding to the mass of the 


object producing the field. 
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lewis clearstrom’(.) that — and aH are Killing vectors 
of g (generating time translations and ¢- rotations respectively). 
Denote by E ("energy at infinity") and L ("angular momentum") the 


corresponding conserved quantities, and let 


BE AS 
MG M 
(2) 
cad ola le 
a aay 
Orbital motion takes place in a plane; choose this plane to 
be 
T 
(3) 6 = 5 


The equations of the orbit are then 


ala m= a 


(4) (Gy? = 8 - Per) 
with 

or eNO 
(5) Hr) = 1a-—9a+E? , 

1s 
(6) => 
r 
and 
| dt E 

(7) at : 2r 

De 


V(r) defined by (5) is an "effective potential" for the orbital motion 


of the test body. The curve parameter tT is proper time for the body. 
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Stable circular orbits occur at the minimum of Vand are 


“xa 


only possible for L > = ri for example 
3 


6r 
(8) ——— 
Gea a Pe Yale 


gives the value of r for a body moving with angular momentum L ina 
circular orbit. Taking r= R in (6) and (7) determines > and t as 


linear functions of T. 


We shall therefore confine our attention to the orbit 
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the orbit is also described by ie To Calo ay =O 
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From (12) we find (cf. 5:.4(7/)) 
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The tangent field to (9) is 
(14) tT > (v"(t)) = (k,0,0,2) 
The nonvanishing Christoffel symbols are 
0 -t -2 Q iL rs —-2 -a 
T's = ror e ; lO Set etete 
1 eae! US | eee 
M4 tort ae A P59 re 
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G15) P33 = =revs sineo , lh5 =r 
ee = -sin @ cos @ , ae = eb 
Jee 
P53 = cot 6 
where 
oe 2x, 
(16) a ee 


has been introduced for brevity. Next, the nonvanishing coordinate com- 


ponents of the Riemann curvature tensor are: 
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However, of more importance are the curvature components 
a q in a Lorentz frame (e,)- First of all, in the natural static 
c 


orthonormal frame 


oie 





the curvature components are: 


(21) 


i 

5 
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A general frame (e,) is obtained from (e) by a Lorentz transfor- 
e 


mation (L*b): 


/ &22) e =e ie and 0° = ie 9° 
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where 


Then for components in the general frame (e,); 


(23) Reet eR ee oe ee 
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(b) Basis for approximation. In terms of any suitable norm l| + | 


we see from (21) that 
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The extension of L to L®L@®LO@L acting on the right on the space 


of (1,3)-tensors has norm 


oy i Q 2 
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The assumption is now made that for r >> Lop ||L|| is of the order 1 
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(meaning that the body is not moving too fast with respect to the 
source r= 0). Note that all considerations are made and, strictly 


speaking, valid here only in Schwarzschild coordinates. 


To obtain a criterion for approximation, we first re-examine 
AG G21 is 
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O 2 pvys ~o 
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v6 
RaBys © 
iteisenowsgeneral| yeconcededs ((12)551131,0127),)129))"that-in 
a curved space-time momentum, p” = Mee » need not be parallel to the 
velocity vo , although this is more of a formal distinction than prac- 
tical ([13]). In practical situations, the approximation is made where- 


by terms containing quadratic spin factors are neglected. However, we 


base our approximation on the assumption that 
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(28) r 


Using (27) we have, independent of the frame ew 


mM? s"r 
(29) A = aus tae (ce=" le 
r 


In the usual weak-field approximation, r >> tO and therefore (28) is 
satisfied provided the spin is not too large. (Note that the various 
discussions ([14], [4], [43]) of the satellite experiment assume “A his 


negligible.) 
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For the earth-gyroscope system mentioned earlier, assume that 


the gyroscope is a spherical 100 gm. mass with a 2 cm. radius and an 


angular velocity of wW = 10° ras 7 and that ates sorbitine 500 kilometers 


above the earth. We temporarily revert to units in which the speed of 


light has its usual value: c= 3x107° cm Hae . Then 


s+ 1.6x10° em ant See 
— See 13 
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M~me= 100 gm R~ 10" cm 


Ca eee =o oct ans 
so that 
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Thus, for a reasonable physical situation A is extremely small com- 


pared to unity. 


(c) Approximate solution for The stirst thing to note as 
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that if the procedure used to solve 84.6(10) is applied to the system 


§5.6(13) and §5.6(14) we obtain 
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where A is given in (27). This shows that both uv and ve are of 
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Equations, 5.o(26) to 5.8(3) to Lowest order are: 
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Using these approximate values in 5.8(23) and 5.8(24) we get 
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Finally, substitution of (34) into 5.8(16) yields after some manipula- 


tion 
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(Compare this with 5.7(23).) 


We shall use (35) to give an estimate of the spin precession; 


LES ty 


dS d aA A a a _A B 
— = ay 
eo 2) Sale deans pea 4 Z 3 dQ ; 


whence 


ee A BLA a Oh La os 
5 = e549 ) oN B en ag en I B io) 5 


so that using the values 5.8(7, 9, 11), 


oO O Se CR aeNe} Nee @ 
(36) (Se + 252g,)8 = e, (8 ve + ku 


Substitution of 5.8(8,10) together with (35) into (36) gives 


oh a *B Mieke fo} ene 2,-2.0 
= + + 
(37) (5. SE Sea (Ss Boe Me By6 en MS “S Va) 
Bo,72, ~0 -4_2~ 
x § Meee S Ge eM oe 
( oe o 
where 
v =iwt sR p 
or} a LVvVea 
POPS Woe oe elk 0.0 
(38) Ceeere Gr oto V 
[hh el yPy? 
OL 2 uvea|o 
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so that the spin precession is second order in i , and therefore 


undetectable in any practical sense. 


We now estimate the extent to which S% undergoes Fermi- 
Walker transport along the orbit (recall 84.6e). The definition 
4.6(45) of i: is compatible with the definition (28) for  (c = 1) 
in orders of magnitude. Comparing (39) with 4.6(31) (they differ by 
Oe equation 4.6(52) shows that Fermi-Walker transport of spin is 
first order in \ for the free case as well for constrained motion. 
Equation (31) says that for constrained motion, Fermi-Walker transport 


is given by 4.6(52) corrected by second-order terms. 


~fRb ¥ 








409°") fee , © of tee ee es c 


ABH Ub “Sq oa RMIAD wi! ovonkseg iene aw ¥ 
nolsinPué ed? .{sa- bh Tie shies ads ea dee to 

= 5) * x98 (8S) sotzte2ien air Mw olacsegee are * ila 

| a21fRb yoke) (oeye. 4 ieee OBE) ambaemats shu gtingats 3. 680 tk 
iiqa 30 toneeshis teolieW—bbiesF veut ewiiia (seas s ooliswps cone 

tzom legtlatdeans yO). Tee Re Sao eee oe wa) A Be seein a 
THqe a7 | Ji<tiew? ood Ion beeda eau = a exes (Ie). netsaypl a 
“190 ‘gthegubgiebin at hetenten (S#)0.4 4d cove at sea 


Le 


* 


> 






iy 7. 


1) 


2) 


3) 


4) 


=) 


6) 


7) 


8) 


ues 


Footnotes 


Often we interpret a vectorfield & e¢ ¥(M) as a linear differ- 


ential operator of c (M), the ring of all C -functions on M. 


In this discussion "e" no longer denotes the identity of the 


group. 


A symplectic form is a closed 2-form which is non-degenerate at 
each point; a presymplectic form is a 1-form whose exterior 


derivative is a symplectic form. 


Equation (33) shows that for free fall motion in curved space 
Euler's equations remain the same as in flat space. On the other 
hand, geodesic motion of the centre of mass is destroyed by the 


curvature of 2%. We meet an equation like (34) in 84.6. 


This is a notational device only. In particular, a term including 
v as a factor need not make sense in the Galilei case - it is 


merely ignored. 


The more suggestive term "integrable" is beginning to replace 
&& 3 & 


Ufiret, order lat’. 


One often sees the definition that a contact structure is a 1-form 


age 


-8@ such that 6 A (a8)* is a volume on M where r = ~—(dimM —- 1). 


2 


In our context, 98 is called an exact contact structure ([1]). 


But more general is Souriau's model of a spinning particle in 
Minkowski space ([48], page 181). The evolution space is diffeo- 
4 Zi 


MOEPHLG tO Ki <x Re x 5? Whichetsenot —[*0 for any Q;5 “and & is 


not exact. 
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Chosen to be future-pointing with respect to the orientation. 


Other than finite-dimensional. 


Not to be confused with the term as used in the theory of ordinary 


differential equations. 


A worldline is a 1-dimensional timelike submanifold of V. 
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GLOSSARY OF NOTATIONS 


The following is a list of commonly used symbols and their 


meanings as used in the thesis. Reference to the text is given where 


appropriate. 

General. 

= definition 

siti Tf andvonly 1c: as equivalent. to 
E ApeonvermAs i secontainedsins) 

€ ae A, a is an element of A 
Maps. 


few AP Beye a > f(a) .peemaps each “a <A to f(a) eB 


£ (A) image of A under f 

ao inverse image of S cB: Pee) = {ave A= f(a) € St 
fog composition, g followed by f 

= inverse map to f , ote = id 


Linear Algebra. 


Vows ess finite-dimensional vector spaces 
R real numbers 
Ro = RX, 6 XR » n times; prototype n-dimensional vector space 
a , : a 
(A » matrix whose (a,b)-entry is A b 
a b ab 


Note: AP , A ee Ab ; A, eA weal denote an (a,b) entry 
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lea! frame of V 
a 
components; i EV 
k 
Vv dual space to V 
# 
{97} frame of V dual to e, 
a 
OF Kronecker delta 
ne transpose of A 
# 
Rog natural pairing between V and V 
(ri) an inner product on V 
® tensor product 
A exterior product 
4 X +4 w, interior product of the p-form w by the vector X , 


yielding a (p-1) - form 


Lie Theory. 

Gal aL tes eroupsm(o2.se) 

(Sigg Ne ee their Lie algebras 

(eee bracket operation in g 

GL(n) general linear group of nxn invertible real matrices (§2.3g) 
g2(n) Lie algebra of GL(n) 

{EP} standard basis of g(n) 

SO (3) orthogonal 3X3 matrices with positive determinant (82.4a) 

so (3) Lievalgebra of 50 (3) eU32.4c) 

Gal homogeneous Galilei group (84.2a) 


Lor connected component of identity of Lorentz group (84.2b) 
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Manifolds. 

ee continuously differentiable to any order; smooth 

Juha Ny ee manifolds: ol , Hausdorff, paracompact, connected (§2.3a) 
nbd neighbourhood 

(U,h) local chart at meM determining local coordinates (x) 


k 
TMs TM tangent space, cotangent space at meM (32 33b) 


#* 
Tas ol aM Gee Me? eM tangent and cotangent bundles of M 
Eve ce Cc (M) real-valued smooth functions on M 
Xe *(M) smooth vector fields on M _ (§82.3a) 
* 
o ¢ ¥ (M) smooth 1-forms on M 
W € #? (M) smooth p-forms on M 
{o, 2 ae local C (M) - module basis of ¥(M) 
Ox 
a co * 
{dx } local C (M) - module basis of #X# (M) 


Maps and Operations on Manifolds. 


ie Ee re smooth mapping between manifolds 


f£. : IM? TN induced differential map (82.3b) 


* 


x Ox * 
£ +TN-23TM induced dual map (§2.3b) 


* 
is » ¥P (nN) +#? (mM) pullback of p-forms 


me : . : 
clan: x? (Mm) > x? 1 oy) exterior derivative of forms 
Le Lie derivative by the vector field xX 
covariant derivative in the 0. - direction 
oe Qa 
ov covariant derivative along a curve 
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Frame Bundles (also see the Appendix). 


L(M) bundle of linear frames of M 

SO (M) bundle of orthogonal frames (83.2b) 

tory) .Galc) bundles of {Lor |Gal} frames on space-time V (84.3b) 
wo g - valued connection 1-form on G(V) 

A fundamental (Killing) vector field associated to Ae g 


(xv ,e0) coordinates on L(M) , Lor(V) , or Gal(V) 


{e.sEnoEa} global basis of either ¥(Gal(V)) or #(Lor(V)) 
Miscellaneous. 

ibe a Lagrangian (§2.5b;83.3a) 

V 3-dimensional gradient, Vf = (0,f) 

m (rest) mass of a body 

By ihe rotational and translational kinetic energy (83.3a) 
M mass Or apspinning body. Mea ami Vv eres (84.5d) 

n Minkowski metric, n = diag(-1,1,1,1) 

g Lorentz metric on space-time V 

EABC 3-dimensional alternating tensor 

Eupys Levi-Civita density 

(ae) symmetrization of indices 


[ ] antisymmetrization of indices 
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APPENDIX 


Principal Bundles, Frame Bundles and Connections 


(a) Group actions. Let G be a Lie group and M amanifold. A 


left actionfof 8Guton "™ Visa Ctha map 


(1) @ ?,GxXM 29M. (a,x) a2 ax 


Wy 


Y(a,x) 


such that (ba)x = b(ax) for all a,beG and x€ M. A right action 


is similarly defined as J) : MXG>M. Sometimes the notation 
RO) = W(x,a) = xa 
(2) 
Meee dows. =) = be 
= 
is used. 


We say that G acts effectively on M LEE RG) =x for 
all xeM implies a=id. G acts freely on M iff there is an 


x € M such that RG) =x implies a = id. 


(b) Principal fibre bundle. Let M be a manifold and G a Lie 
group. A principal fibre bundle over M with group G is a manifold 


P. together with an action of BGiten F satisfying: 
(i) G acts freely on P on the vight: 


Re PkCo Se (usa) = R,(u) = ua 
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(ii) M is the quotient space of P by the equivalence relation 
induced by G , and the canonical projection T7 : PUM eas 


smooth. 


(417) Ser Sis locally trivial; that is, about each x ¢€M there is a 
neighbourhood U_ such that ro is isomorphic to UxXG in 
the sense that there is a diffeomorphism jp : a +> UxG 
such that w(u) = (T(u),¢(u)) , where > : 7 (WU) + G satis- 


fies (ua) = (6(u))a for all ue n+) And ale Go. 


Cali PP" the bundle space,” i the base space, G the struc— 
ture group, T7 the projection and ae the fibre over x = M-. 


~ 


(Note that ra) = (Ce formaliuex co Mi) 


If M is any manifold and G a Lie group, then trivially we 
havessee=iIMxG -with the action) PXG > FP , (Gan Dee ab) asean 


example of a principal fibre bundle. 


(c) Fundamental vector fields. Let P be a principal bundle with 
group G and base M. The action W of G on P induces an isomor- 
phism o of the Lie algebra ¢ of G@ nto (Pe) Li) Ave ¢7 then 


exp(tA) is a 1-parameter Subgroup of) Ge, at ue Fo define 
d 
(3) OO ede tony gee. 


* 
o(A) is variously denoted in the literature by A Cie alge: 


(Sometimes A is called a "Killing vector field".) 
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* 
The mapping A>A 9 open get Sue TP is an isomorphism of g 


Wathe thestangent space-at\u to the fibre through uw : ~g ~ Tm (m(u)) 


(d)""* Connection on a principal fibre bundle. “Let P be a principal 
bone Lesover=aMoe With croup meGe.m Lig Wo cube 5 CG, denotes the subspace of 
TP tangent to the fibre a rey - Ap connection 9) in P is ’a smooth 
assignment of a subspace Q0 of the tangent space to each u€ P_ which 


satisfies: 
(i) TP = CG. @ QF 


EY Way ate © 


a u 


Gy is called the vertical subspace and OF the horizontal subspace of 
TP Te AnVeector in Gy is called vertical, while one in oF isecaived 


horizontal. 


(e) Connection form. Let I be a connection on P and denote 
by v(X) (resp. h(X)) the vertical (resp. horizontal) component of X 


with respect to I. Define the g - valued connection 1-form w on P 


as follows: if xX « TO then 


w(X) =A where A_ satisfies 
(4) 


ne 2 y@ 


A is necessarily unique (since g = GC. JP andi OOO) e= 20 eth Ses 


horizontal. 
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(f) Bundle of linear frames of a manifold. Let M be a n - dimen- 
sional manifold. A linear frame at x e¢M is an ordered basis 
e = (Ejo+++5e) of TM . Denote by LM the set of all such frames of 


TM -) Then UCM) = v LM is a principal fibre bundle over M with 
xeM 


DEOTectIOn il: (*,e) > x= and a rightsaction ote clin) bys) ((x,e)5 2). > 


e e = 3 
(x,e°a) where (e a), age ic 


Let (x) be a local coordinate system on a neighbourhood U 


Che xe Mo. Then each frame ~@ of TM can be uniquely expressed as 
(5) Revere noe (e-) e'Gi (a) 
i il 3c : a 


‘ a a 
We take, then, as a local coordinate system at (x0 eye Cue (x »e,) on 


the neighbourhood n+ (U) 


Occasionally, we view a frame e at x asa non-singular map 


Re TM . The action of GL(n) on L(M) is then the composition 


ie Bias ye —S TM for ae GL(n) and ee LM : 


(g) The canonical RnR” - valued 1-form § on L(M) is defined by 
a(x) = e (my (X)) 
(6) COe=Fce Gin 


CO er Et, cle = (x, ec) ands ee is interpreted as in the last para- 
u 


F oh ue! . 
graph. In a local coordinate system (x 7c) » we can write: 


See m4) 
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(h) A connection [ on L(M) is called a linear connection on M . 


Let {E} be the standard basis of g(n) ; then in terms of a local 


coordinate system (x"',e0) on L(M) , the g&(n) - valued connection 1- 


form w has the expression 


(8) 


a b 
W = W b Eo 
a _ pa } a Y B 
Wy 664 a, + a By e. dx” ) 
a depends on x = 71(x,e) only. 


(i) Standard horizontal vector fields. Given a linear connection 


on M anda vector €€ Ro ,» define the horizontal vector field € by 


(9) 


Seon 


T(E] ,) = e(&) 


: youl = Kx, eye (Meee: Roe TM 


& is called the standard horizontal vector field corresponding to 
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Lae ie In a local coordinate system (x"',e,) : 
VS 0! re) Go oO 
Se So eae ve Sa ae 
dx Bee 
a 
The standard horizontal vector fields corresponding to the 
é n joe ee) 
standard basis (S.) of 8 ((S.) = 6) ; 
) op) 9 ae) 
(10) ez=e.(—-M ee) , 
ig Yr ee QO a de. 


forms a basis of the horizontal subspace QF for each u = (x,e) € L(M) 


(j) Parallelism of ThM) . Recablithat tto each Ae g2(n) there 


is associated the fundamental vector field A which spans the vertical 


subspace G, at each ueéL(M) . Locally, we have 
a Gd 
oe ap “b Lo 
de 
a 


Thus, the aia vectors {e | E| } forms a basis of 
Mali bp) u 


T LOD foreach | ulce Liha wt) Xee T Lao , then 


a a —b 
(Ga) bE + UL EA 
, ; Onlin au 
determines local coordinates (x »e05V »U,) for TL(M). ~ Furthermore, 


the following orthonormality conditions hold: 


I 
Oo 


a Bo Dae 
a” (e,) a(E.) = 0 


(13) 


Cc 


a a Cees Ane 
®.(e.) = 0 Ww pee = 84 5, 


bese 


ghiell twice Lathes! 25a ishawse OT 7 , is Gi 7 
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(k) Structure equations of a linear connection. Define a 2-form 


QOeeon LG) by: 


(14) 2) (XX) = dw Ch (OK) nh (X)) 


: ; =a. Db 
where X,>X, e T LO) and h() is the horizontal part. (= os EA 


is called the curvature 2-form of the connection and satisfies the struc-— 


ture equation 


(15) ae ee a ie Ga 


The canonical 1-form § (equation (7)) satisfies a second 


structure equation: 
(16) de> + uw A ge = O° 


where O(X, »X,) = d6(h(X,) ,h(X,)) 4g the torsion form of [2 in 


this thesis, we assume Gr=vOos. (But see Trautman [Sa\ons | S61)- 
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